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ABSTRACT 


In  its  motion  response  to  a  disturbance,  a  ship  can  be  closely  ap- 
proximated as  a  second  order  system.  The  solution  of  these  second  order 
equations  is  used  to  determine  the  motion  stability  of  the  ship.  By  ap- 
plying a  step  disturbance  to  a  ship  model  and  measuring  the  response, 
the  linearity  can  be  used  to  find  the  response  to  any  disturbance. 

In  this  work  a  step  disturbance  was  applied  to  a  Mariner  class 
hullform  and  the  response  measured.  From  this  response  the  response  to 
regular  sinusoidal  excitations  was  computed.  The  response  to  the  sinu- 
soidal inputs  was  then  used  to  determine  the  coefficients  of  virtual 
inertia  and  damping  for  the  Mariner  in  pitch  and  heave. 

The  derived  results  at  zero  speed  were  compared  with  values  com- 
puted by  strip  theory, 

THESIS  SUPERVISOR:  J,  E.  Kerwin 

TITLE:   Associate  Professor  of  Naval  Architecture 


-l- 


ACKNOWLEDGEMENTS 

The  authors  wish  to  thank  Professor  J.  E,  Kerwin  for  his  continual 
advice  and  encouragement  during  this  project. 

The  authors  also  wish  to  thank  Mr.  Frank  Sellars  for  his  support, 
especially  in  obtaining  the  analog  to  digital  conversion  of  the  data. 

The  assistance  of  Commander  W,  R,  Porter  was  highly  beneficial  in 
obtaining  the  theoretical  results. 

The  assistance  of  Mr.  H,  der  Kinderen  in  performing  the  experi- 
mental work  was  greatly  appreciated. 

This  work  was  done  in  part  by  the  computation  center  at  the  Mas- 
sachusetts Institute  of  Technology,  Cambridge,  Massachusetts, 


-11- 


TABLE  OF  CONTENTS 

Page 

ABSTRACT  i 

ACKNOWLEDGEMENTS  ii 

TABLE  OF  CONTENTS  iii 

INDEX  OF  TABLES  v 

INDEX  OF  FIGURES            ¥  vi 

CHAPTER  I  INTRODUCTION  1 

CHAPTER  II  PROCEDURE  5 

2.1  Model  Characteristics  5 

2.2  Test  Conditions  5 

2.3  Step  Disturbance  Application  6 
2.U  Measurement  of  the  Step  Response  11 
2.5  Test  Procedure  12 

CHAPTER  III  ANALYSIS  lli 

3.1  Background  lh 

3.2  Static  Coefficients  15 

3.3  Digitizing  Rate  17 
3.ii     Theoretical  Analysis  19 

CHAPTER  IV     RESULTS  23 

ii.l     Heave-No  Pitch  Condition  23 

U.2     Pitch-No  Heave  Condition  25 

Iu3     Theoretical  Results  26 

CHAPTER  V     DISCUSSION  OF  RESULTS  52 

5.1  Experimental  Errors  52 

5.2  Character  of  the  Results  5k 

5.3  Theoretical  Results  56 

-iii- 


Page 

CHAPTER  VI  CONCLUSIONS  £8 

6.1  Evaluation  of  Step  Method  £8 

6.2  Experimental  Results  59 

6.3  Theoretical  Results  59 
CHAPTER  VH  RECOMMENDATIONS  61 

7.1  Experimental  Technique  6l 

7.2  Future  Studies  63 
BIBLIOGRAPHY  6k 
APPENDIX  *  A-l 
APPENDIX  A  A-2 
APPENDIX  B  A-l* 
APPENDIX  C  A-6 
APPENDIX  D  A-8 


3  IV- 


INDEX  OF  TABLES 


Table  I  Characteristics  of  the  Model 

Table  II  Static  Restoring  Coefficients  in  Heave 

Table  HI  Effect  of  Convolution  Interval  on  Coefficients 

Table  IV  Section  Transformation  Coefficients 

Table  V  Virtual  Mass  Coefficient  in  Heave 

Table  VI  Damping  Coeff icieiit  in  Heave 

Table  VH  Virtual  Inertia  Coefficient  in  Pitch 

Table  VHI  Damping  Coefficient  in  Pitch 

Table  IX  Theoretical  Virtual  Mass  and  Damping 

Table  X  Heave  Rod  Friction 

Table  XI  Transformation  Coefficients 

Table  XII  Comparison  of  Errors  with  Those  Found  by 

Process  Used  in  the  Main  Body  of  the  Project 


Page 

16 

19 
21 

27 
30 
33 
36 

39 

A-3 

A-ll 

A-12 


-v- 


INDEX  OF  FIGURES 

Page 

Figure  I  Pitch-No  Heave  Configuration  7 

Figure  U  Electric  Wiring  Arrangement  9 

Figure  HI         Magnet  Parting  Circuit  9 

Figure  IV  Heave-No  Pitch  Configuration  10 

Figure  V  Virtual  Mass  in>  Heave  as  a  Function  of 

Froude  Number  1±0 

Figure  VI  Virtual  Mass  in  Heave  at  Zero  Froude  Number  lil 

Figure  VTI        Damping  in  Heave  as  a  Function  of  Froude 

Number  hZ 

Figure  Viil   Damping  in  Heave  at  Zero  Froude  Number  U3 

Figure  IX     Virtual  Inertia  in  Pitch  as  a  Function  of 

Froude  Number  hk 

Figure  X     Virtual  Inertia  in  Pitch  at  Zero  Froude 

Number  U5 

Figure  XI  Damping  in  Pitch  as  a  Function  of  Froude 

Number  ii6 

Figure  XII    Damping  in  Pitch  at  Zero  Froude  Number  hi 

Figure  XHT   Comparison  of  Theoretical  and  Experimental 

Zero  Speed  Virtual  Mass  in  Heave  U8 

Figure  XIV    Comparison  of  Theoretical  and  Experimental 

Zero  Speed  Damping  in  Heave  h9 

Figure  XV    Comparison  of  Theoretical  and  Experimental 

Zero  Speed  Virtual  Inertia  in  Pitch  50 

Figure  XVI    Comparison  of  Theoretical  and  Experimental 

Zero  Speed  Damping  in  Pitch  £l 

Figure  XVII   Zero  Speed  Heave  Response  (Parallel)  A=5> 

Figure  XVTII  Zero  Speed  Heave  Response  (Perpendicular)  A-£ 

Figure  XIX    Comparison  of  the  Actual  Mariner  Station  One 

with  the  Calculated  Section  A-7 

-vi- 


CHAPTER  I 
INTRODUCTION 

The  equations  of  ship  motion  developed  by  Abkowitz  (1)  are  a 
linearized  set  of  second  order  partial  differential  equations.     In  the 
vertical  plane  the  equations  are 

(m-Z..)z    -  Z-z   -  Z  z    -    (Zv+mX_)e--    (Z.+u  Z..)9   -    (Z  +u  Z.  )6  =  Z(ex)  (l) 

v        z '  z  z  VS       G  0     o    z  xQoz 

(I  -M..)6   -    (M-+U  M   )9   -    (M  +u  M.  )0  -    (M,.+nK_)z"  -  M.z   -  M  z   =  M(ex)        (2) 
v   y      g'  x    Q     o   zy  v    0     o   z'  v   z        G  z  z 

z  =  displacement  in  the  vertical  plane 

z  s  velocity  in  the  vertical  plane 

*z  *  acceleration  in  the  vertical  plane 

9  =  angular  displacement  in  pitch 

0*  ■  angular  velocity  in  pitch 

8  ■  angular  acceleration  in  pitch 

u0  ■  forward  velocity  of  the  ship 

m  ■  mass  of  the  ship 

Z"    =  added  mass  coefficient 

z 

Z-  ■  damping  coefficient 

Z   =  static  restoring  force  coefficient 

X.,    B     distance  from  origin  to  the  center  of  gravity( equal  to  zero) 

G 

Z-  "  coupled  virtual  mass 

8 

Z*  =  rotary  force  derivative  coefficient 

Z   ■  static  force  derivative  coefficient 

8 

I   =  moment  of  inertia 

y 

M--  a     added  moment  of  inertia  coefficient 
9 

M«  ■  damping  coefficient 
8 
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Mfi    ■     static  restoring  moment  coefficient 

M-    ■»     coupled  virtual  moment  of  inertia  coefficient 

M^     «     rotary  moment  derivative  coefficient 

M      ■     static  moment  derivative  coefficient 

z 

These  equations  show  a  coupling  effect  between  pitch  and  heave. 
However,  by  decoupling  these  equations  and  neglecting  second  order 
effects,  these  equations  can  be  reduced  to 

(m  -  Z.-)z  -  Z.z  -  Z  z  =  Z(ex)  (3) 


(Iy  -  M~)e  -  M^6  -  MQ9  =  M(ex)  (k) 


The  equations  are  now  in  a  form  that  is  readily  solvable.  The 
solution  of  these  uncoupled  equations  is  highly  significant  in  deter- 
mining the  motion  characteristics  of  a  ship.  These  characteristics  are 
important  in  determining  stability  of  motion  as  well  as  factors  affect- 
ing the  comfort  and  safety  of  personnel  embarked  on  the  ship.  Before 
these  equations  can  be  solved,  however,  it  is  first  necessary  to  deter- 
mine the  hydrodynamic  coefficients  of  added  inertia  and  damping. 

The  added  inertia  is  a  term  that  results  from  the  fact  that  a  ship 
moves  in  a  dense  medium  and  therefore  must  displace  a  heavy  mass  of 
liquid  whenever  it  moves  in  any  plane.  This  movement  of  the  mass  of 
liquid  has  the  effect  of  essentially  increasing  the  mass  of  the  ship. 
The  damping  is  a  term  associated  with  energy  dissipation  and  not  vis- 
cous effects.  As  a  surface  ship  operates  at  the  interface  between  water 
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and  air,  any  vertical  motion  results  in  the  creation  of  waves.  The 
energy  lost  in  wave  creation  acts  as  a  damping  force  to  reduce  the  os- 
cillations of  the  ship. 

Many  methods,  both  theoretical  and  experimental,  have  been  de- 
veloped to  determine  these  coefficients.  One  such  method  was  developed 
by  Narita  (11)  at  the  Massachusetts  Institute  of  Technology  in  1°62. 
It  consisted  of  applying  a  known  moment  step  disturbance  to  a  hull  form 
and  measuring  the  transient  response  in  pitch,  while  restraining  the 
model  from  motion  in  heave.  The  resultant  motion  of  the  model  was  a 
transient  wave  of  a  decaying  sinusoid  form  leading  to  the  step.  When 
this  transient  response  was  obtained,  it  was  used  to  compute  the  theo- 
retical  response  of  the  model  to  any  known  input.  This  makes  use  of  a 
property  of  a  linear  system,  that  when  its  response  to  a  step  or  im- 
pulse input  is  known,  its  response  to  any  other  input  may  be  found  by 
use  of  the  convolution  integral.  As  a  ship  is  closely  approximated  by  a 
linear  second  order  system,  this  principle  is  applicable. 

This  property  can  be  used  to  apply  known  regular  sinusoidal  inputs 
and  determine  the  ship's  response.  This  in  essence  produces  the  same 
results  as  could  be  obtained  by  making  forced  oscillation  tests.  The 
output  of  the  regular  sinusoidal  inputs  can  then  be  used  to  calculate 
added  mass  and  damping. 

Narita  (11)  computed  added  inertia  and  damping  for  the  pitch 
motion  only.  In  this  work  the  procedure  was  extended  to  measure  the  co- 
efficients for  pitch  and  heave.  Two  series  of  tests  were  run5  one  in 
which  the  model  could  pitch  with  heave  restrained  and  one  in  which  the 
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model  could  heave  with  pitch  restrained. 

The  authors  know  of  no  previous  tests  conducted  on  the  Mariner 
hull  form  for  the  purpose  of  determining  coefficients  of  damping  and 
added  mass.     However,  a  theoretical  prediction  for  the  zero-speed  con- 
dition may  be  approximated.     Porter  (lU)   has  reported  that  the  conformal 

transformation  of  the  unit  circle  represented  by  the  equation 

N 

.-(2n+l) 


n=0 

where  z  is  the  plane  of  the  actual  ship  section  and  £  is  the  plane  of 
the  unit  circle,  can  represent  a  given  station  of  the  ship.  N  is  se- 
lected and  the  coefficients  a2n*i  are  found  which  give  the  best  fit  to 
the  actual  section.  Lewis  (10)  used  this  transform  with  N  ■  1  and 
Landweber  and  Macagno  (9)  extended  it  to  N  *  2.  In  the  present  work, 
a  method  was  used  which  can  give  a  least  mean  square  error  approximation 
to  the  actual  ship  section  if  the  number  of  coefficients  is  selected. 
The  least  number  of  coefficients  which  gave  a  reasonable  fit  for  each 
section  was  used. 

After  the  coefficients  of  the  transform  were  determined,  Porter's 
method  was  used  to  obtain  coefficients  of  damping  and  added  mass  for 
each  station  at  several  frequencies  and  the  results  were  numerically  in- 
tegrated throughout  the  length  of  the  model  by  using  strip  theory  to 
obtain  added  mass  and  damping  for  the  model. 
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CHAPTER  II 
PROCEDURE 

2.1  Model  Characteristics 

The  tests  were  conducted  using  a  model  of  the  Mariner  class  hull 
form,  which  has  the  characteristics  shown  in  Table  I. 

table  1^ 
CHARACTERISTICS  OF  THE  MODEL 


Length  between  perpendiculars 

Beam 

Draft 

Displacement  in  fresh  water 

Block  coefficient 

Radius  of  gyration 

Static  restoring  moment  (pitch) 

Static  restoring  force  (heave) 

Material  of  model 

2.2  Test  Conditions 


5.28  feet 
.76  feet 
.27  feet 
h2.h   pounds 
.6125 
1.26  feet 
l£8. Qk   foot  -  pounds/ radian 
II18.8O  pounds/foot 
Fiberglass,  reinforced  plastic 


The  tests  were  conducted  at  the  Massachusetts  Institute  of  Tech- 
nology ship  model  towing  tank.  The  tests  were  divided  into  two  series, 
pitch  with  no  heave  and  heave  with  no  pitch.  Each  series  consisted  of 
tests  conducted  at  Froude  numbers  of  0.0,  0,1,  0.2,  0.3.  Further,  the 
tests  conducted  at  zero  Froude  number  were  conducted  parallel  to  and 
perpendicular  to  the  centerline  of  the  towing  tank  in  order  to  investigate 
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the  effects  of  wall  reflections.  In  the  pitch  -  no  heave,  configuration, 
the  ballast  weights  were  shifted  in  order  to  determine  the  effect  of 
changing  the  radius  of  gyration  on  the  computed  coefficients. 

2.3  Step  Disturbance  Application 
2.3*1  General 

As  the  tests  were  to  be  conducted  in  a  liquid  medium,  it  was  de- 
cided that  the  high  potential  devices,  as  recommended  by  Narita  (11), 
posed  a  danger  to  personnel  working  with  the  equipment  and  to  other 
electronic  equipment  in  the  event  of  a  wiring  mistake  or  a  faulty  lead. 
For  this  reason  low  current  and  voltage  devices  were  investigated,  and 
it  was  decided  that  the  best  device  might  be  a  small  electromagnet.  The 
magnet  chosen  for  these  tests  was  a  simple  iron  core  wound  magnet.  For 
the  purposes  of  these  tests,  the  magnet  required  less  than  two  amperes 
of  current  and  less  than  two  volts  of  potential.  These  values  were  not 
only  safe  to  personnel  in  the  event  of  a  casualty,  but  were  easily  ob- 
tainable with  standard  direct  current  power  supplies.  Initially  bat- 
teries were  employed,  but  as  the  tests  took  a  long  time,  because  of  de- 
lays in  waiting  for  the  tank  to  settle  to  a  dead  calm,  the  batteries 
rapidly  lost  their  potential.  In  place  of  the  batteries  a  Sorenson 
Model  QR  18-1. 5>A  direct  current  power  supply  was  used. 

2.3,2  Pitch -No  Heave  Configuration 

To  apply  the  pitch  disturbance,  the  model  was  mounted  as  shown  in 
Figure  I.  A  small  metal  plate  was  mounted  on  the  bow  and  on  the  stern 
of  the  model.  A  millimeter  bar  was  used  to  mount  the  magnet.  This 
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Figure  I 
PITCH-NO  HEAVE  CONFIGURATION 


Towing  Carriage 


Millimeter  Bar 
Electromagnet 


n. 


Heave  Rod 
Pitch  Bearing 
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enabled  the  operators  to  have  greater  ease  in  setting  the  disturbance  at 
very  close  to  one  and  a  half  degrees.  As  mentioned  above,  the  current 
was  supplied  to  the  magnet  from  a  Sorenson  Model  QR  18-1, 5A  direct  cur- 
rent power  supply.  The  electric  wiring  arrangements  were  as  shown  in 
Figure  II.  By  releasing  the  current  to  the  magnet,  the  ship  will  let  go 
under  a  step-like  displacement,  following  a  decaying  sinusoid  coming  to 
rest  at  a  zero  angle  of  pitch. 

As  the  exact  moment  is  desired  at  which  the  ship  actually  reacts 
to  the  disturbance,  a  second  circuit  was  constructed  as  shown  in 
Figure  HI,  When  the  ship  begins  to  react,  the  magnet  parts  from  the 
steel  plate  and  produces  a  step  reproduction  on  the  recorder.  The  be- 
ginning point  of  this  step  can  be  used  to  determine  the  starting  point 
of  the  model's  response  and  thus  increase  the  accuracy  of  calculation 
based  on  the  response  of  the  model. 

2.3.3  Heave  -  No  Pitch  Configuration 

In  constructing  the  heave  arrangements,  it  was  necessary  to  devise 
a  method  of  lifting  the  entire  ship.  It  was  also  desirable  to  not  in- 
troduce any  forces  tending  to  cause  pitch,  which  would  thus  put  undue 
strain  on  the  heave  rod.  For  this  reason  it  was  decided  to  mount  a 
small  metal  plug  directly  on  the  heave  rod  and  to  lift  the  model  by 
connecting  the  magnet  to  the  metal  plug.  As  the  heave  rod  is  mounted 
approximately  at  the  longitudinal  center  of  floatation  of  the  model,  es- 
sentially no  pitch  moment  is  introduced.  The  magnet  was  fastened  to  the 
carriage  by  two  aluminum  rods  as  shown  in  Figure  IV.  The  circuitry  is 
exactly  the  same  as  for  the  pitch  configuration. 
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Figure  II 

ELECTRIC  WIRING  ARRANGEMENT 
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MAGNET  PARTING  CIRCUIT 
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Figure  IV 
HEAVE-NO  PITCH  CONFIGURATION 
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In  both  heave  and  pitch  testing  the  magnet  and  metal  connecting 
plate  were  mounted  rigidly,  so  there  was  no  concern  for  catching  a  fall- 
ing weight.  Then  to  reconnect,  it  was  only  necessary  to  close  the 
switch  and  lift  the  model  so  that  the  plate  and  the  magnet  came  in  con- 
tact, at  which  time  the  magnet  would  support  the  weight  of  the  model. 
Then  it  was  only  necessary  to  wait  until  the  tank  settled  to  an  absolute 
calm  before  beginning  the  next  run, 

-a 

2,U  Measurement  of  the  Step  Response 
2,U,1  Pitch 

In  order  to  increase  the  validity  of  linearity  of  motion  theory, 
an  initial  displacement  of  only  one  and  a  half  degrees  was  used.  To  get 
the  maximum  accuracy  in  measuring  this  small  variation,  two  methods  were 
investigated,  use  of  a  linearsyn  and  use  of  a  pitch  bearing.  To  use  a 
linearsyn,  maximum  excursion  of  motion  is  desired.  This  would  be  at  the 
ends  of  the  model.  However,  the  motion  here  is  circular  enough  to  cause 
non-linear  motion  of  the  linearsyn  and  create  inaccuracies.  Also,  the 
vibrations  near  the  ends  of  the  model  were  found  to  be  the  greatest,  and 
this  too  would  have  an  adverse  effect  on  measurements  made  by  a  linear- 
syn, 

A  pitch  bearing  located  right  at  the  longitudinal  center  of 
floatation  has  the  advantage  that  it  can  measure  angular  motion  right 
at  the  point  about  which  pitch  takes  place.  The  pitch  bearing  is  a 
rotary  motion  displacement  pickoff  of  a  differential  transformer  type 
and  is  highly  sensitive  to  any  rotary  displacement.  Therefore,  it  was 
decided  to  employ  the  pitch  bearing  for  the  pitch  measurements, 
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2,Iu2  Heave 

To  measure  heave,  the  linearsyn  was  perfectly  suited.  It  is  a 
linear  motion  displacement  pickoff  of  the  differential  transformer  type. 
It  consists  of  a  coil  assembly  and  movable  magnetic  core  and  is  capable 
of  detecting  motions  as  small  as  0.000001  of  an  inch  between  the  core  and 
coil  assembly.  Mounted  directly  onto  the  heave  rod  assembly,  the 
linearsyn  gives  an  accurate  reading  as  to  the  movement  of  the  model  in 
the  heave  plane. 

To  increase  the  validity  of  model  linearity  in  heave  motion,  a 
disturbance  of  only  four- tenths  of  an  inch  was  applied  to  the  model. 
This  proved  sufficient  to  obtain  several  well-defined  cycles  of  oscil- 
lation* 

2,5  Test  Procedure 

After  the  model  was  mounted  on  the  carriage  and  circuitry  arranged, 
the  magnet  was  connected  to  either  displace  the  model  in  pitch  or  heave, 
depending  on  which  tests  were  being  conducted.  Tests  were  then  carried 
out  in  two  groups,  pitch  and  heave,  at  Froude  Numbers  of  0.0,  0,1,  0,2, 
0.3. 

With  the  model  displaced  it  was  necessary  to  wait  until  the  tank 

i 
had  settled  to  an  absolute  calm  before  beginning  a  test.  Then  the  model 

was  brought  up  to  speed  and  the  switch  opened,  stopping  power  to  the 

magnet.  At  that  instant  the  ship  would  fall  toward  zero  in  a  decaying 

sinusoidal  motion,  characteristic  of  a  second  order  system.  This  motion, 

as  well  as  the  exact  instant  at  which  the  magnet  parted  from  the  metal 

plate,  was  recorded  on  both  the  Sanborn  recording  paper  and  on  magnetic 
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tape  for  later  analysis.  To  begin  the  next  run,  the  model  had  only  to 
be  towed  back  to  the  opposite  end  of  the  tank,  the  magnet  re-engaged 
with  the  metal  plate,  and  then  allowed  to  remain  suspended  until  the 
tank  settled. 

For  the  zero  speed  runs  the  model  was  placed  at  the  center  of  the 
tank,  first  parallel  to  and  then  perpendicular  to  the  centerline.  Here 
the  procedure  was  the  same  as  before  except  that  there  was  no  forward 
motion  of  the  model. 
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CHAPTER  III 
ANALYSIS 

3.1  Background 

The  analysis  procedure  used  In  this  work  was  essentially  that  used 
by  Kerwin  and  Narita  (£)•  A  step  input  is  applied  to  a  system,  and  re- 
sponse is  represented  by 

P(t)  -i(l  +R(t))  (6) 

where  C  is  either  the  static  restoring  force  in  heave  or  the  static  re- 
storing moment  in  pitch,  R(t)  represents  the  transient  portion  of  the 
response  and  l/C  is  the  final  value  the  response  seeks  after  the  tran- 
sient has  died  out. 

To  get  the  response  to  a  regular  sinusoidal  input,  the  convolution 
integral  is  used. 

V0(t)  -  J    V  (t)  H(t  -  t)  dT  (7) 

H(t)  represents  the  impulse  response  of  the  system.  However,  as  the  step 
response,  S(t),  is  known,  it  may  be  used  with  the  input  in  the  form 

r°°         dV.  (T) 

v*>-L  -v-s(t-T)dT  (8) 

and  the  same  results  will  be  obtained.  From  this,  the  response  of  the 
model  to  sinusoidal  inputs  in  either  heave  or  pitch  may  be  obtained. 

To  obtain  the  damping  and  virtual  mass  coefficients,  we  represent 
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our  second  order  equations  as  follows % 


*r  J  +  bp  i  +  cp6=  Mo  slnoJt  (9) 

d.  z        dz 

^H  77  +     \  dt   +  CHZ=  Zo   Slnajb  (10) 

at 


The  solution  of  these  equations  will  be  of  the  form 

a   s incut  +  b   cosyjt  (ll) 

from  which 

.  Cp-Vpf/^/)  (12) 

(JU 

-M  b 

w(ap    +  bp  ) 

.     CH   -  VA2   ^H2)  (14) 

-Z  b 

O^H      +  bH    ) 


3«2  Static  Coefficients 

Equation  (6)  indicates  that  the  step  response  tends  to  l/c  as  the 
transient  dies  out.  In  the  initial  simplification  of  the  equations s   it 
was  assumed  that  this  was  a  static  term  and  not  speed  dependent,  Ger- 
ritsma  (2)  shows  that  this  term  is  speed  dependent  and  that  this  depend= 
ency  is  further  a  function  of  the  hull  form.  While  the  variation  in 
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speed  is  not  great  at  very  low  Froude  numbers,  it  does  approach  approxi- 
mately 10  percent  at  a  Froude  number  of  three-tenths,,  This  is  still  not  a 
large  effect,  but  it  may  exert  a  measurable  influence  on  the  values  of 
added  mass  and  damping  calculated  at  the  various  Froude  numbers. 

Accordingly,  tests  were  conducted  to  determine  the  change  in  the 
static  restoring  coefficients  with  speed.  Tn  pitch  a  known  moment  was 
applied  to  the  ship  and  the  pitch  angle  was  measured.  As  the  pitch 
velocity  and  acceleration  were  *both  zero,  the  static  restoring  coeffi- 
cient could  be  found  by  dividing  the  moment  by  the  resultant  angle.  The 
static  restoring  moment  at  zero  speed  was  computed  to  be  15>8.8U  ft. -lb./ 
radian.  This  same  procedure  was  repeated  at  Froude  numbers  of  0.1,  0.2, 
and  0.3,  and  a  variation  of  less  than  1  percent  was  observed  at  each 
speed.  Therefore,  the  static  coefficient  at  zero  speed  was  used  for  the 
analysis  at  each  of  the  Froude  numbers. 

A  known  heaving  force  was  then  applied  to  the  model,  and  heave  was 
measured  at  Froude  numbers  of  0.0,  0.1,  0.2,  0.3.  In  this  case  the  re- 
sults followed  the  form  of  those  found  by  Gerritsma  as  shown  in  Table  U„ 


Table 

II 

STATIC 

RESTORING  COEFFICIENTS 

IN  HEAVE 

Froude  Number 

Zz 

0 

li;8.8  lbs/ft 

.1 

lla.O  lbs/ft 

.2 

1U2.8  lbs/ft 

.3 

139.2  lbs/ft 
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The  variation  is  not  large,  but  is  significant  enough  to  warrant 
using  the  value  found  at  its  corresponding  Froude  number, 

3,3  Digitizing  Rate 

The  data  as  obtained  from  the  experiments  was  recorded  in  analog 
form  on  magnetic  tape.  As  the  convolution  integral  was  performed  by 
series  methods,  it  was  necessary  to  have  the  data  in  digital  form.  Pre- 
vious experiments  conducted  by  Kerwin  and  Narita  (5)  showed  the  results 

-a 

of  added  mass  and  damping  to  be  highly  dependent  on~the  digitizing  rate 
used,  with  damping  being  the  more  critical.  From  their  tests  of  a 
known  second  order  system,  it  was  found  that  a  digitizing  rate  of  one 
kilocycle  per  second  gave  exact  results  for  added  mass  and  to  within 
three-tenths  of  1  per  cent  for  damping. 

Accordingly,  a  digitizing  rate  of  one- thousandth  of  a  second  was 
used.  To  digitize  the  data  at  this  rate,  it  was  necessary  to  record  a 
one  kilocycle  per  second  square  wave  on  the  tape,  as  the  process  used 
with  the  IBM  709U  computer  reads  the  value  of  data  everytime  there  is  a 
positive  to  negative  vertical  change  in  the  digitizing  signal. 

As  mentioned  earlier,  the  response  of  the  model  was  recorded  as 
well  as  a  signal  that  showed  the  instant  the  magnet  parted.  These  sig- 
nals were  recorded  first,  and  the  square  wave  was  recorded  afterward. 
The  square  wave  was  recorded  insuring  that  it  started  prior  to  the  in- 
stant of  the  magnet  parting  signal.  By  this  procedure  and  obtaining  the 
digitized  data  for  both  the  response  data  channel  and  the  parting  signal, 
the  parting  signal  data  could  be  scanned,  and  the  exact  instant  of  mag= 
net  parting  could  be  found.  On  the  other  channel,  this  was  the  starting 
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point  for  which  the  response  signal  data  was  used  in  the  analysis  pro- 
gram. 

As  mentioned  above,  to  get  maximum  accuracy,  a  digitizing  rate  of 
a  thousand  points  per  second  was  used  to  obtain  the  best  representation 
of  the  model  response  data.  In  convolving  this  data  with  the  input 
sinusoid,  a  series  method  was  used.  This  involved  shifting  the  sinu- 
soid in  reference  to  the  response  data.  The  response  to  the  sinusoid 
was  computed  as  a  series,  and  the  final  value  summed.  In  shifting  the 
waveform,  it  was  found  that  the  fine  spacing  of  .001  seconds  was  not 
required.  The  analysis  program  was  altered  to  determine  what  interval 
could  be  used  and  still  produce  the  same  results  that  one-thousandth  of 
a  second  produced.  The  main  purpose  of  this  alteration  was  to  speed 
up  the  analysis  program  and  reduce  computer  time  used  in  obtaining  re- 
sults. Table  IH  shows  the  results  of  increasing  the  convolution 
interval  at  several  frequencies  or  corresponding  periods  of  excitation. 

These  values  were  obtained  using  ideal  data  for  the  second  order 

2.5  ^#  +  7.5  §  +  1776  (16) 

dt 

From  Table  HI  it  can  be  seen  that  increasing  the  convolution  in- 
terval actually  improved  the  accuracy  of  results  at  all  except  the  high- 
est frequency,  period  equal  to  0.6U  seconds.     In  this  case  the  added  in- 
ertia stayed  approximately  the  same,  but  the  damping  became  worse  by 
seven-thousandths.     The  over-all  results  did  show  that  larger  convolution 
intervals  gave  accurate  results.     Accordingly,  larger  intervals  were  used 
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in  the  analysis,  depending  on  the  input  excitation  frequency. 


Table  in 

EFFECT  OF  CONVOLUTION  INTERVAL 

ON  COEFFICIENTS 

Period 
(Seconds) 

Convolution  Interval 
(Seconds) 

Added  Inertia 
(Slug  Feet*) 

Damping 
(Slue  Feet 
Second 

1.28 

.001 

2.50.7 

7.581 

1.28 

.010 

2.503 

7.550 

1.28 

.010 

2.503 

7.5U9 

1.28 

.020 

2.503 

7.51i9 

1.28 

.032 

2.502 

7.5U6 

1.28 

.oho 

2.502 

7.51i5 

2.56 

.001 

2.511 

7.658 

2.56 

.010 

2.506 

7.602 

2.56 

.016 

2.506 

7.601 

2.56 

.020 

2.506 

7.600 

.61* 

.001 

2.U98 

7.529 

,6k 

.010 

2.U99 

7.536 

.6U 

.020 

2.1*99 

7.536 

3.U     Theo: 

retical  Analysis 

In  the  present  problem  offsets  for  21  stations  of  the  Mariner  hull 
were  read  off  a  ^-inch-scale  body  plan,  and  by  iteration  on  the  IBM  709U 
computer,  the  coefficients  a2n+i  which  gave  the  best  fit  to  the  actual 
sections  were  determined.  No  more  than  five  such  coefficients  for  any 
one  station  were  used  because  the  physical  significance  of  coefficients 
beyond  &g  becomes  unclear.  The  values  of  the  coefficients  giving  the 
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best  fit  for  each  section  are  tabulated  in  Table  IV. 

After  the  coefficients  *2n+l  were  determined,  they  were  used  as 
inputs  to  a  computer  program  which  utilizes  Porters  method  as  explained 
in  reference  lii  and  calculates  coefficients  of  damping  and  added  mass 
for  each  station.  The  results  for  each  station  were  integrated  over  the 
length  of  the  model  by  using  strip  theory  to  give  the  coefficients  for 
the  model. 
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Table  TV 
SECTION  TRANSFORMATION  COEFFICIENTS 


Station 

Half -Be am 
0.20 

* 

Draft 
27.0 

Section  Area 
Coefficient 

0 

7.1407 

1 

U.30 

27.0 

0.999 

2 

9.75 

27.0 

0.818 

3 

16.15 

27.0 

0.779 

h 

22.70 

27.0 

0.791 

5 

28.50 

27.0 

0.830 

6 

33.25 

27.0 

0.861* 

7 

36.30 

27.0 

0.911* 

8 

37.70 

27.0 

0.958 

9 

38.00 

27.0 

0.982 

10 

38.00 

27.0 

0,988 

11 

38.00 

27.0 

0.979 

12 

38.00 

27.0 

0.965 

13 

38.00 

27.0 

0.932 

U* 

37.75 

27.0 

O0868 

15 

36.1*5 

27.0 

0.787 

16 

33.1*0 

27.0 

0.700 

17 

28.25 

27.0 

0.609 

18 

21.05 

27.0 

0.519 

19 

12.20 

27.0 

0.390 

20 

1.80 

0.85 

0.020 
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Table  IV  (continued) 


Statior 

\            al 

a-3 

!£ 

!Z 

!£ 

0 

-.85755187 

-.11659060 

-.0189161*91* 

+.0065186719 

-.00037253161* 

1 

-.66039k 

-.102210 

-.006052 

+.02111*2 

2 

-.1+7391+1 

-.01*8562 

+.012796 

+  .01721*6 

+.0061+57 

3 

-•273985 

+  .021*895 

+  .021*931 

♦.015372 

•■•• 

1* 

-•112^73 

-.017126 

+.019732 

+  .0071+28 

+.007161 

5 

♦.001188 

-.032628 

+.019083 

+.000588 

+.005890 

6 

+.076083 

-.0531+92 

+.016873 

+.000628 

+  .005291+ 

7 

+.112537 

-.078719 

+.010882 

-.001901 

+.001656 

8 

+  .11+3558 

-.103009 

+  .001*781* 

9 

+  .1501+10 

-.118791 

-.001681 

-.002359 

10 

+.150708 

-.125100 

-.003798 

-.002210 

♦ .000777 

11 

+.151138 

-.123101 

-.003217 

-.003000 

12 

+.15081*3 

-.110398 

-.001079 

-.001*638 

13 

+  .11*8500 

-.085690 

+  .0051*38 

-.001*678 

ll+ 

+  .11*0856 

-.01*81*39 

+.015873 

-.001*178 

+.000558 

IS 

+  .121*071* 

~. 002855 

+  .0231+61 

-.001276 

♦.000786 

16 

+  .0811*37 

+  .01*9610 

♦.029596 

-.001735 

17 

-.022515 

+.101750 

+  .01*2262 

-.0001*02 

+.005171 

18 

-.220882 

♦.138063 

+.078202 

+.001977 

+.001509 

19 

-.551*295 

+.IIO698 

+.101173 

+  .035001* 

+.020560 

20 

+.319277 

-.081*215 

+  .051+789 

-.052018 

*. 019595 
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CHAPTER  IV 
RESULTS 

i 

ii.l  Heave -No  Pitch  Condition 

To  obtain  values  of  added  mass  and  damping  in  heave,  four  tests 
were  made  at  each  of  the  Froude  numbers  desired.  A  typical  response 
curve  is  shown  in  Appendix  B  ,  These  curves  are  presented  to  show 
wall  reflection  effects,  but  they  also  show  the  typical  response  at  zero 
speed.  At  forward  speeds  greater  than  zero  the  model  response  curves 
were  much  the  same,  but  considerable  vibration  was  also  introduced  be- 
cause of  the  towing  carriage  oscillations. 

li.1.1  Virtual  Mass 

The  tests  conducted  at  each  speed  were  averaged  to  obtain  the  vir- 
tual mass.  This  average  value  was  then  nondimensionalized  by  dividing 
it  by  the  mass  of  the  model.  Each  of  these  values  is  shown  in  Table  V 
as  a  function  of  excitation  frequency.  Furthermore,  each  of  the  ob- 
served values  was  compared  with  the  average  to  obtain  the  worst  excur- 
sion. This  worst  excursion  from  the  average  is  expressed  as  a  percent 
of  the  average  in  Table  V. 

In  most  cases  all  points  were  vexy  close  together,  and  percentage 
errors  were  small.  However,  occasionally  a  derived  value  was  at  vari- 
ance with  the  rest  of  the  results  and  its  error  appeared  quite  large. 
This  took  place  at  those  values  farthest  from  the  model1 s  natural  fre- 
quency of  approximately  7.6  radians  per  second. 
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The  symbols  used  in  Table  V  are  as  follows; 
w  radian  frequency  of  excitation 
A  virtual  mass 
A1  nondimensional  added  mass 
g  gravity 
A  displacement  of  the  model 

The  average  values  of  the  virtual  mass  are  plotted  in  Figure  V 
for  each  of  the  four  Froude  numbers.  Figure  VI  is  a  comparison  of  the 
zero  speed  values  obtained  with  the  model  parallel  to  and  perpendicular 
to  the  centerline  of  the  tank* 

iul.2  Damping 

The  same  procedure  was  used  with  damping  as  was  used  with  added 
mass.  The  results  of  the  average  values,  as  well  as  the  maximum  percent 
error,  are  shown  in  Table  VI.   The  damping  values  were  nondimensional- 
ized  by  multiplying  the  damping  obtained  by  the  square  root  of  the 
length  of  the  model  times  gravity  divided  by  the  displacement  of  the 
model. 

Errors  in  damping  were  very  small,  but  were  generally  slightly 
higher  than  for  added  mass*  However,  near  the  natural  frequency  of  the 
model,  the  errors  were  very  small,  increasing  as  the  frequency  was  in- 
creased or  decreased. 

Additional  symbols  introduced  in  Table  VI  areg 
B  damping 

B'  nondimensionalized  damping 
L  length  of  model 
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The  average  values  of  damping  are  plotted  in  Figure  VTI  for  each  of 
the  four  Froude  number's.  Figure  VIII  is  a  comparison  of  the  zero  speed 
values  obtained  with  the  model  parallel  to  and  perpendicular  to  the  cen- 

V 

teiiine  of  the  tank. 

1|.2     Pitch-Mo  Heave  Condition 

To  obtain  the  values  of  the  added  inertia  and  damping  coefficients 
in  pitch,   six  tests  were  conducted  at  each  of  the  Froude  numbers,,     Three 
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tests  were  conducted  with  the  bow  initially  up,  and  three  tests  were  con- 
ducted with  the  stern  initially  up.  The  response  curves  were  the  same  as 
for  the  case  of  heave,  in  that  each  was  a  decaying  sinusoid.  At  speeds 
other  than  zero  there  was  a  large  amount  of  vibration  in  the  response^ 
paralleling  the  situation  found  in  heave. 

U.2.1  Virtual  Inertia 

The  same  analysis  program  was  used  to  compute  the  value  of  virtual 
inertia  in  pitch  as  was  used  to  compute  the  virtual  mass  in  heave.  The 
only  changes  were  to  change  several  constants  and  the  nondimensionalizing 
factor.  Virtual  inertia  was  nondimensionalized  by  multiplying  the  value 
found  by  the  value  of  gravity  divided  by  the  displacement  times  the 
length  squared.  The  values  of  virtual  inertia  found  at  each  of  the  four 
Froude  numbers  are  shown  in  Table  VTI,  In  this  table  A  designates  the 
value  of  virtual  inertia.  Figure  XX  is  a  plot  of  the  virtual  inertia 
versus  excitation  frequency  for  each  of  the  four  Froude  numbers.  Fig- 
ure X  is  a  comparison  of  the  zero  speed  condition  with  the  model  paral- 
lel to  and  perpendicular  to  the  centerline  of  the  tank.  This  figure 
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also  shows  the  effect  on  the  virtual  inertia  of  changing  the  radius  of 
gyration  of  the  model  from  1,26  feet  to  1.17  feet. 

li.2.2  Damping 

The  damping  was  computed  using  the  same  analysis  program  as  for 
heave  with  the  appropriate  constants  changed.     It  was  nondimensionalized 
by  multiplying  by  the  square  root  of  the  length  times  the  value  of  grav- 
ity and  dividing  this  square  root  by  the  displacement  times  the  length 
squared.     The  values  of  damping  as  a  function  of  Froude  number  and  fre- 
quency are  shown  in  Table  VTTI.     The  damping  is  plotted  in  Figure  XI  to 
show  its  variation  with  Froude  number.     Figure  XH  is  a  plot  of  the  zero 
speed  conditions,  also  showing  the  effect  of  the  reduced  radius  of  gyra- 
tion. 

U.3     Theoretical  Results 

Based  on  the  coefficients  of  transformation  as  presented  in  Table  IV, 
Porter's  theory  was  applied  to  determine  the  coefficients  of  damping  and 
virtual  mass  in  both  pitch  and  heave.     These  coefficients  were  determined 
at  nine  frequencies  for  each  station,  and  the  results  were  numerically 
integrated  by  strip  theory  to  obtain  the  coefficients  of  virtual  mass  and 
damping  for  the  model.     These  results  are  presented  in  Table  IX.     These 
same  results  are  plotted  in  Figures  XHI-XVI  to  compare  their  values  with 
the  experimentally  obtained  zero  speed  values. 
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Table  V 


VIRTUAL  MASS  COEFFICIENT  IN  HEAVE 


(1)     Froude  Number 

0.3 

w 

A 

U.13k 

1.951 

h.909 

1.7l4l 

6.01*2 

1.821* 

7.11*0 

1.788 

7.851* 

1.752 

8.727 

1.718 

9,817 

1.702 

11.220 

1.759 

12.083 

1.81*0 

13.090 

1.891 

(2)     Froude  Number      0.2 


w 

A 

A' 

U.131* 

3.01*5 

2.313 

I*.  909 

2.308 

1.751 

6.01*2 

2.129 

1.617 

7.11*0 

1.959 

1.1*86 

7.851* 

1.838 

1.395 

8.727 

1.732 

1.316 

9.817 

1.651 

1.253 

11.220 

1.601 

1.217 

12.083 

1.625 

1.231* 

13.090 

1.678 

1.273 

A1  Max.  %  Error 

1.1*81  10.08 

1.322  8.73 

1.386  1.01* 

1.358  2.13 

1.330  1.1*8 

1.303  .93 

1.292  .65 

1.336  1.1*2 

1.398  2.1*1* 

1.1*37  3.75 


Max.  %  Error 
11*.  86 
9.50 
8.65 
6.1*8 
5.33 
1*.27 
1*.18 
3.26 
7.08 
8.11 
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Table  V  (continued) 


(3)     Frcude  Number 

0.1 

w 

A 

A1 

lul3U 

3.U76 

2.619 

U.909 

2.692 

2.0ia 

6.0U2 

2.282 

1.73k 

7.11l0 

2.113 

1.605 

7.85U 

1.978 

1.501 

8.727 

1.815      . 

1.1*01 

9.817 

1.757 

1.332 

11.220 

1.71k 

1.301 

12.083 

1.75k 

1.331 

13.090 

1.862 

l.iak 

(10     Froude  Number 

0.0  (Parallel) 

w 

A 

A« 

k.l3k 

3.876 

2.938 

ii.909 

2.8k5 

2.159 

6.01*2 

2.273 

1.723 

7.II4O 

2.186 

1.659 

7.85U 

2.05k 

1.555 

8.727 

1.923 

1.U59 

9.817 

1.870 

I.I4I8 

11.220 

1.905 

l.hkh 

12.083 

1.879 

l.k22 

13.090 

1.993 

1.511 

Max.  %  Error 
li.21 
2.71 
7.50 
5.35 
2.68 

3.63 
U.10 
5.5U 
5.70 
3.kk 

Max.  %  Error 

.10 
1.23 
2.86 
2.06 
1.80 
1.35 
1.07 

.68 

.30 

.71 
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Table  V  (continued) 


(5)     Froude  Number 

0.0  (Perpendicular) 

w 

A 

A' 

k.llk 

2.923 

2.220 

U.909 

2,696     • 

2.0U2 

6.0li2 

2.267 

1.719 

7.1U0 

2.097 

-» 

1.590 

7.85U 

2.030 

1.539 

8,727 

1.958  * 

1.U8U 

9.817 

1.900 

l.liiiO 

11.220 

1.901 

i.Ua 

12.083 

1.9k9 

1.U78 

13.090 

2.077 

1.572 

Max.  %  Error 
6.88 
7.03 
U.10 

3.39 
ii.28 

U.75 
5.21 
6.81* 
6.57 
3.90 
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Table  VI 
DAMPING  COEFFICIENT  IN  HEAVE 


(1)     Froude  Number 

.3 

w 

B 

B1 

U.13U 

7.333 

2.253 

U.909 

6,180 

1.900 

6.01*2 

5.991 

i.8ia 

7.1U0 

6.k9h 

1.99k 

7.85U 

6.637 

2.039 

8.727 

6.655 

2.01*5 

9*817 

6,1*20 

1.973 

11.220 

5.929 

1.821 

12.083 

6.173 

1.897 

13.090 

7.295 

2.21*1 

(2)  Froude  Number   .2 


w 

B 

B« 

hnUh 

1.8U2 

1.1*88 

k.909 

5.279 

1.621 

6.01*2 

lu910 

1.509 

7.11*0 

5.913 

1.828 

7.851* 

6.213 

1.910 

8.727 

6.281 

1.931 

9.817 

6.182 

1.901 

11.220 

5.636 

1.731 

12.083 

5.138 

1.578 

13.090 

U.871* 

1.U99 

Max.  £  Error 
3.67 
1.38 
1*.62 
2.37 
3.92 
3.61 
3.12 
li.00 
11.12 
9.27 

Max.  %  Error 
10.78 
15.1*8 
15.32 
11.88 

7.62 

1*.00 

8.27 
16.20 

19.70 
100.00 
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Table  VI  (continued) 


(3)     Froude  Number 

a 

w 

B 

B« 

k.ljk 

5.105 

1.661; 

Ii.909 

7.302 

2.21*5 

6.01*2 

6.I29 

1.976 

7.UiO 

6.832 

2.100 

7.85U 

7.063 

2.172 

8.727 

7.0U5 

2.163 

9.817 

6.U63 

1.987 

11.220 

5.2W 

1.612 

12.083 

fc.305 

1.322 

13.090 

1.1*93 

1.380 

(h)     Froude  Number 

0  (Parallel) 

w 

B 

B« 

U.13U 

5.911 

1.812 

U.909 

8.109 

2.582 

6.0U2 

6.660 

2.01*3 

7.lli0 

7.126 

2.186 

7.85U 

7.390 

2.267 

8.727 

7w03U 

2.159 

9.817 

6.112 

1.876 

11.220 

5.550 

1.702 

12.083 

h.991 

1.531 

13.090 

3.U31 

1.052 

Max,  %  Error 
17.67 
lii.00 

7.11 

5.11* 

3.88 

1.11; 

1.89 

9.91 
16.00 
27.21 


Max.  %  Error 
6.96 

5.51 
2.12 

.1*1* 
.03 
.111 

.1*9 

.80 

.20 

U.20 
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Table  VI  ( continued) 


(5)  Froude  Number 

0  (Perpendicular) 

w 

B 

B« 

h.l3U 

5,325 

1.635 

h.909 

7.U30 

2.281 

6.0i*2 

8.1*61 

2.600 

7*11*0 

8,30U 

2,551 

7.85U 

8,2l*2 

2.530 

8.727 

8.005 

2,1*58 

9.817 

7.279 

2.236 

11,220 

5.806 

1.782 

12,083 

h.113 

1.10*8 

13.090 

3,773 

1.158 

Max.  $  Error 
10,13 
10.52 

7,97 

1.13 

1.72 

1.72 
10,86 
22,10 
26.52 
U5.30 
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Table  VII 


VIRTUAL  INERTIA  COEFFICIENT  IN  PITCH 


(1)     Froude  Number       .3 


w 

A 

iul3lt 

2.1*10 

h.909 

2,501 

6.0U2 

2.286 

7.110 

2.538 

7.85U 

2.521* 

8.727 

2.503 

9.817 

2.622 

11,220 

2.399 

12.083 

2.593 

13.090 

2.1*20 

(2)     Froude  Number 

il 

w 

A 

U.13U 

2.986 

U.909 

2.528 

6.0U2 

2.598 

7.1U0 

2.667 

7.85U 

2.605 

8.727 

2.558 

9.817 

2.593 

11,220 

2.lbl* 

12.083 

2.55U 

13.090 

2.U03 

A' 
.066 
.068 
.062 
.069 
.069 
.068 
.071 
.065 
.071 
,066 


A' 
.081 
.069 
.071 
.073 
.071 
.070 
.071 
.066 
.069 
.065 


Max.  $  Error 
7.31; 
6.U0 

3.15 
1.03 
1.38 
1.68 
3.38 
3.71 
U.32 
6.32 


Max.  %  Error 
2.11 
1.U2 
2.38 
1.35 
1.11 

.78 

.10 

.70 
ii.67 
5.86 
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Table  VII  (continued) 

(3)     Froude  Number       .1 

w                                         A  A'                               Max.  %  Error 

U.13U  U.020  .109  h.(& 

h.909  3.538  .096  3.7k 

6.0U2  3.019  .082  U.U5 

7.1U0  2.715  .07U  3.17 

7.85k  2,551  .069  2.70 

8.727  2.1a2     -  .066  U.02 

9.817  2.316  .063  6,09 

11.220  2.200  .060  7.31 

12.083  2.150  .058  9.35 

13.090  2.181  .059  9.81 


A 

A' 

1+.020 

.109 

3.538 

.096 

3.019 

.082 

2.715 

.07U 

2,551 

.069 

2.1412        r 

.066 

2.316 

.063 

2.200 

.060 

2.150 

.058 

2.181 

.059 

[  (Parallel) 

A 

A' 

3.623 

.099 

3.591 

.098 

3.167 

.086 

2,821 

.077 

2.656 

.072 

2.519 

.068 

2.I4I8 

.066 

2.279 

.062 

2,199 

.060 

2.1*0 

.058 

w                                         A  A1                                Max.  %  Error 

U.13U  3.623  .099  .la 

ii.909  3.591  .098  #i|2 

6.01+2  3.167  .086  .36 

7.11+0  2.821  .077  .32 

7.85U  2.656  .072  .31 

8.727  2.519  .068  .21 

9.817  2.ia8  .066  .17 

11.220  2.279  .062  .32 

12.083  2,199  .060  .10 

13.090  2.11+9  .058  .19 
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Table  VII  (continued) 


(5)     Froude  Number 

0  (Perpendicular) 

w 

A 

A' 

ii.l3i| 

1.999 

.05U 

U.909 

3.296 

.090 

6.0U2 

3.017 

.082 

7.1U0 

2.622 

.071 

7.85k 

2,U80 

.067 

8,727 

2.1*08     - 

.065 

9.817 

2.U91 

.068 

11.220 

2.273 

.062 

12.083 

2.3h5 

.06U 

13.090 

2.553 

.069 

(6)     Froude  Number 

0  (Reduced  R) 

w 

A 

A» 

U.13U 

3.079 

•  08U 

1.909 

3.062 

.083 

6.0U2 

2.67U 

.073 

7.11*0 

2.357 

.06U 

7.85U 

2.235 

.061 

8.727 

2.13U 

.058 

9.817 

2.07k 

,056 

11,220 

2.037 

,055 

12.083 

1,998 

.05U 

13.090 

2.0ii9 

.056 

Max.  £  Error 
3,10 
1.U0 
1.28 
1.37 
1.59 
1.7U 
2.69 
lu53 
3.67 
1.61 


Max.  %_  Error 

.29 

,ii6 

.10 

.1*7 

.76 
1.03 
1.30 
1.67 
2.20 
2.30 
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Table  VIII 


DAMPING  COEFFICIENT  IN  PITCH 


(1)     Froude  Number 

.3 

w 

B 

U.13U 

6.651 

h.909 

8.5H* 

6.Ql*2 

6,157 

7.11*0 

5.887 

i.m 

5.876 

8.727 

5.155 

9.817 

5.130 

11.220 

l*.3l*3 

12.083 

2.021 

13.090 

It.  762 

(2)     Froude  Number 

si 

w 

B 

l*.13l* 

8.521* 

U.909 

7.821* 

6.01*2 

5.1*68 

7.11*0 

5.863 

7.8^1* 

5.823 

8.727 

5.189 

9.817 

li.981 

11.220 

3.588 

12.083 

2.921 

13.090 

2.003 

B'  Max.  %  Error 

.073  32.30 

.09U  11,30 

.068  13.95 

.065  8.92 

.065  li.95 

,057  2.26 

.057  2.50 

.01*8  52.60 

.022  100.00 

.052  iiO.50 


B'  Max.  #  Error 

.091*  1*.1*2 

.086  2.58 

.060  6,70 

.065  3.H* 

.061*  .86 

.057  1.00 

.055  3.1*8 

.01*3  13.60 

.032  21.60 

.022  100.00 


-36- 


Table  YIU  (continued) 


(3)      Froude  Number 

0.1 

w 

B 

B» 

Max.  %  Error 

U.13U 

3.087 

.03U 

1U.95 

ii.909 

U.316 

.0U8 

1.01 

6.0li2 

5.35h 

.059 

1.36 

T.lbO 

6.0U8 

.067 

2.1*8 

7.851* 

6.318 

.070 

1.2k 

8.727 

5.900      , 

.065 

3.72 

9.817 

5.U51* 

.060 

5.81 

11.220 

h.665 

.051 

9.60 

12.083 

3.328 

.037 

13.36 

13.090 

1.276 

.oil* 

20.82 

(h)     Froude  Number 

0  (Parallel) 

w 

B 

B1 

Max,  %  Erroi 

l*.13l 

.820 

.009 

.51 

h.909 

3.99k 

.ohk 

1.0U 

6.0U2 

5.508 

.061 

.1*2 

7.11*0 

6.173 

.068 

.53 

7.85U 

6.091 

.067 

.38 

8.727 

5.637 

.062 

.55 

9.817 

5.102 

.056 

2,66 

11,220 

U.205 

.0li6 

2.72 

12.083 

3.155 

.035 

2.U0 

13.090 

1.U91 

.016 

ll*.l*0 
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Table  VHI  (continued) 


(5)     Froude  Number 

0  (Perpendicular) 

w 

B 

B1 

l*.13l* 

-2.723 

-.030 

U.909 

1.351 

.015 

6.0l*2 

6.197 

.068 

7.11*0 

7.333 

.081 

7.851* 

7.080 

.078 

8.727 

6.353 

.070 

9.817 

5.552 

.060 

11.220 

7.11*5 

.079 

12.083 

3.531* 

.039 

13.090 

10.385 

.111 

(6)     Fronde  Number 

0  (Reduced  R) 

w 

B 

B» 

l*.13l* 

1.81*7 

.020 

li.909 

3.886 

.01*3 

6.0i*2 

6.085 

.067 

7.1U0 

6.550 

.072 

7.851* 

6.1*55 

.071 

8.727 

6.157 

.068 

9.817 

5.566 

.061 

11.220 

5.395 

.059 

12.083 

lu701 

.052 

13.090 

1*.578 

.050 

Max,  %  Error 

1.1*0 

11*.  36 

2.57 

1,01* 

.10 

.88 

1.15 

9.61* 
10.83 
23.60 


Max,  %  Error 
16.83 
5.90 
2.26 
1.1*3 
1.25 

.99 

.1*1 
1.26 
5.10 
1.75 
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Table  IX 

THEOEETICAL 

VIRTUAL  MASS  AND  DAMPING 

Frequency 

Virtual 

Mass 

Damping 

Virtual 
Inertia 

Damping 

1**11 

2.823 

8.U31 

U.OlU 

9.605 

5.82 

2.U2U 

8#5U0 

3.U68 

10.231 

7.13 

2.295 

7.703 

3.230 

9.857 

8,23 

2.257 

<    6.733 

3.139 

9.272 

9.20 

2.260 

5.933 

3.105 

8.536 

10.08 

2.281 

5.050 

3.073 

7.826 

10.89 

2.310 

ii.153 

3.115 

6.378 

11.6U 

2.3ii0 

3.893 

3.13U 

6.588 

12.80 

2.369 

3.086 

3.157 

5.U36 
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FIGURE  V 
VIRTUAL  MASS  IN  HEAVE  AS  A  FUNCTION  OF  FROUDE  NUMBER 
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FIGURE  VI 
VIRTUAL  MASS  IN  HEAVE  AT  ZERO  FROUDE  NUMBER 
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FIGURE  VII 
DAMPING  IN  HEAVE  AS  A  FUNCTION  OF  FROUDE  NUMBER 
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FIGURE  VIII 
DAMPING  IN  HEAVE  AT  ZERO  FROUDE  NUMBER 
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FIGURE  IX 
VIRTUAL  INERTIA  IN  PITCH  AS  A  FUNCTION  OF  FROUDE  NUMBER 
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FIGURE  X 
VIRTUAL  INERTIA  IN  PITCH  AT  ZERO  FH0UDE  NUMBER 
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FIGUHE  XI 
DAMPING  IN*  PITCH  AS  A  FUNCTION  OF  FROUDE  NUMBER 
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FIGURE  III 
DAMPING  IN  PITCH  AT  ZERO  FROUDE  NUMBER 
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FIGURE  XIII 

COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL 
ZERO  SPEED  VIRTUAL  MASTlN  HEAVE 
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FIGURE  XIV 

COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL 
ZERO  SPEED  DAMPING  IN  HEAVE 
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FIGURE  XV 

COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL 
ZERO  SPEED  VIRTUAL  INERTIA  IN  PITCH 
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FIGURE  XVI 


COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL 
ZERO  SPEED  DAMPINGTn  PITCH 
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Chapter  V 
DISCUSSION  OF  RESULTS 

5,1     Experimental  Errors 

The  resultant  curves  of  virtual  mass  and  damping  as  shown  in  Fig- 
ures 7  -  XII  show  the  respective  values  for  the  Mariner,  In  comparing 
these  curves  with  those  of  other  similar  type  large  hullforms,  these 

-a 

curves  show  the  same  trends  in  their  values 0     These  curves  are  gener- 
ally believed  to  give  the  values  for  the  Mariner  at  the  various  speeds 
for  the  given  excitations.     However^,   there  are  several  experimental 
errors  or  limitations  which  may  limit  their  accuracy. 

Initially  the  magnet  was  supported  by  a  wire.     When  the  tests 
were  conducted,   the  random  vibrations  were  so  great  as  to  almost  ob- 
scure the  actual  record  of  the  model's  motion.     For  this  reason  stiff 
supporting  members  were  used.     This  greatly  reduced  the  vibrations 9  but 
the  inherent  vibrations  in  the  towing  carriage  structure  and  members 
could  not  be  eliminated.     These  vibrations  became  greater  at  the  higher 
speeds.     They  were  also  affected  by  the  rate  at  which  the  model  was 
brought  up  to  speedy  being  much  greater  when  the  model  was  quickly  ac- 
celerated to  the  desired  speed.     In  general  the  vibrations  were  high 
frequency  with  little  effect  on  the  results s  but  since  the  vibrations 
were  so  random,   it  is  likely  that  there  were  enough  low  frequency  com- 
ponents to  affect  the  results. 

Another  source  of  error  was  introduced  by  the  procedure  used  in 
recording  the  results.     Initially  the  response  was  recorded  at  3,75 
feet  per  second.     At  this  low  speed  a  one-kilocycle-per°second  square 
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wave  would  not  reproduce  on  the  tape.  In  transferring  the  response  to 
a  tape  at  30  feet  per  second,  it  was  observed  that  noise,  primarily 
sixty  cycles  per  second,  was  picked  up.  The  sixty  cycle  noise  should 
have  had  no  effect  on  the  results,  but  any  lower  frequency  components 
would  have.  This  effect  was  considered  to  be  very  minor. 

As  mentioned  earlier,  a  second  circuit  was  constructed  to  record 
the  actual  instant  at  which  the  model  began  to  drop,  as  evidenced  by 
the  magnet  parting  from  the  metal  plate.  By  digitizing  this  response, 
it  was  hoped  to  be  able  to  pinpoint  the  proper  point  on  the  model  re- 
sponse data  and  not  introduce  any  errors  by  starting  the  analysis  at 
the  wrong  instant.  The  system  was  generally  satisfactory,  but  it  was 
sometimes  difficult  to  be  sure  if  the  chosen  starting  point  was  within 
two  or  three  thousandths  of  a  second  of  the  actual  one.  To  find  the 
effect  of  missing  the  proper  starting  point,  some  of  the  data  was  ana- 
lyzed at  points  away  from  the  initial  point.  The  virtual  mass  coef- 
ficient showed  little  sensitivity  to  small  changes,  while  damping 
tended  to  be  much  more  dependent  on  the  point  at  which  analysis  was  be- 
gun. Thus,  this  is  a  possible  cause  for  the  errors  being  generally 
larger  in  damping  than  in  virtual  mass. 

Because  of  the  use  of  the  rotary  pitch  bearing,  friction  was  neg- 
ligible. In  heave  a  definite  friction  factor  was  present.  However,  as 
discussed  in  Appendix  A  s   the  friction  had  no  effect  on  the  response 
curve  and  thus  virtually  no  effect  on  the  final  values  of  virtual  mass 
and  damping.  To  increase  the  validity  of  the  results,  this  friction 
factor  should  have  been  included  in  this  analysis  program. 
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£.2  Character  of  the  Results 

The  resultant  curves  of  virtual  mass  and  damping  show  the  general 
form  as  expected.  As  predicted  by  Grim,  the  curves  of  virtual  mass 
should  drop  sharply  from  infinity  to  a  minimum  and  then  rise  to  an  as- 
symptotic  value.  The  infinite  point  should  occur  at  approximately  a 

OJU 

radian  frequency  of  — 9.  =  ■£•  .     While  the  range  of  frequencies  investi- 

g 

gated  do  not  show  this  point,  they  do  definitely  show  this  trend  as  it 
is  seen  that  the  curves  tend  toward  infinity  at  frequencies  inversely 
proportional  to  speed,  the  higher  speeds  rising  at  lower  frequencies. 

The  values  of  damping  also  showed  the  same  trend  as  for  other 
ships.  Theory  predicts  that  damping  should  go  to  zero  at  higher  fre- 
quencies. However,  experimental  techniques  have  shown  damping  to  rise 
at  higher  frequencies.  This  indicates  that  the  theory  is  only  good  at 
low  frequencies o  In  some  cases  the  damping  appeared  to  still  be 
falling  to  zero  in  these  experiments.  However,  the  higher  speeds 
showed  a  definite  rise.  The  lower  speeds  were  not  tested  at  higher 
frequencies,  as  this  is  the  general  region  of  towing  carriage  equip- 
ment natural  frequencies.  In  this  region  the  validity  of  the  results 
was  somewhat  questionable  anyway. 

For  the  curves  in  general  a  study  of  the  results  shows  the  agree- 
ment of  the  results  in  the  vicinity  of  the  natural  frequency  of  the 
model  to  be  very  good.  At  frequencies  removed  from  this  value  the  er- 
rors went  up  considerably.  This  indicated  that  this  method  is  accurate 
only  near  the  model's  natural  frequency.  Also,  as  mentioned  earlier,, 
the  errors  in  virtual  mass  were  much  less  than  for  damping,  indicating 
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that  the  procedure  may  be  better  for  determination  of  added  mass  than 
for  damping.  At  the  higher  frequencies  this  effect  became  very  pro- 
nounced but  was  also  affected  by  equipment  natural  frequencies. 

Kerwin  and  Narita  noted  a  sudden  drop  in  damping  at  approximately 
eleven  radians  per  second.  This  was  attributed  to  the  natural  frequency 
of  the  heave  rod.  In  this  series  of  tests  a  much  stiffer  heave  rod  was 
used  in  an  effort  to  extend  the  range  of  frequencies.  This  sharp  dip 
appears  to  have  been  eliminated^  but  there  is  enough  fluctuation  at  fre- 
quencies above  eleven  radians  per  second  to  render  these  results  of 
little  value.  Thus  this  indicates  an  even  stiffer  heave  rod  would  be 
desirable. 

To  overcome  the  problem  of  frequency  range  limitations,  there  is 
also  the  possibility  of  altering  the  model's  natural  frequency.  This 
was  investigated  by  changing  the  position  of  the  ballast  weights  to 
change  the  radius  of  gyration.  In  analyzing  the  effect  of  this  alter- 
ation, the  damping  did  not  appear  to  change  but  fell  somewhat  between 
the  two  other  zero  speed  conditions.  The  virtual  inertia  did  change  as 
expected  and  is  shown  in  Figure  X. 

The  results  obtained  in  pitch  did  not  depend  on  whether  the  model 
was  initially  bow  up  or  stern  up.  In  heave  the  results  were  all  ob- 
tained with  the  model  lifted  up.  If  a  method  were  used  to  put  pressure 
on  the  model  in  order  to  start  with  it  initially  down  in  the  water,  it 
would  probably  produce  the  same  results  as  when  it  started  lifted  up. 
This  supports  the  principle  of  linearity  which  was  assumed  in  this  en- 
tire procedure. 
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£•3     Theoretical  Results 

5.3.1     Transformation  Coefficients 

The  values  of  the  transformation  coefficients  a2n+i  were  found 
which  mapped  a  cylinder  into  the  sections  of  the  Mariner.     The  agree- 
ment between  the  mapped  sections  and  the  actual  sections  was  very  good 
in  the  vicinity  of  amidships,  while  the  greatest  excursions  were  near 
the  ends.     For  example,  at  station  nine  the  maximum  difference  between 
the  mapped  section  and  the  actual  section  was  only  0.08  feet  with  an 
average  excursion  of  0.039  feet.     The  worst  fit  occured  at  station  one, 
which  had  a  maximum  excursion  of  0,9315  feet  and  an  average  excursion  of 
0.ii22  feet.     The  average  excursion  for  the  entire  ship  was  0.168  feet. 
A  plot  showing  the  difference  in  the  mapped  and  actual  form  at  station 
one  is  given  in  Appendix  C. 

3>«3«2     Comparison  of  Theoretical  and  Experimental  Results 

From  Figures  XHI  and  XIV  one  can  see  that  the  theoretical  values  of 
virtual  mass  and  damping  in  heave  compare  quite  favorably  with  the  ex- 
perimentally derived  values.     The  virtual  mass  \ values  show  closer  agree- 
ment than  do  the  damping  values.     However,  the  difference  between  the 
theoretical  and  experimental  damping  is  of  the  same  order  as  in  the  two 
experimental  values  found  from  essentially  identical  response  curves. 
This  shows  the  extreme  sensitivity  of  damping  to  any  minor  differences, 
and  perhaps  strip  theory  is  not  a  precise  enough  method  to  evaluate 
damping,  unless  very  fine  station  spacing  were  used0 

In  pitch  the  agreement  between  the  experimental  and  theoretical 
values  of  virtual  inertia  was  better  than  the  agreement  in  damping. 
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However,   in  both  cases  the  variations  were  greater  than  for  the  corres- 
ponding values  in  heave.    This  variation  was  believed  to  be  the  result  of 
using  strip  theory  on  the  Mariner  hullform.     The  Mariner  has  a  bulbous 
bow  and  a  wide  flared  stern.     These  sections  produced  a  comparatively- 
large  value  for  both  added  mass  and  damping.     To  use  strip  theory  in 
pitch,  each  section  is  multiplied  by  the  distance  from  the  center  of 
gravity  squared  and  then  integrated  over  the  length  of  the  model.     These 
extreme  bow  and  stern  sections  had  a  very  large  arm  to  the  center  of 
gravity  and  thus  produced  very  great  factors.     In  the  actual  case  there 
are  two  reasons  why  the  strip  theory  is  questioned.     Near  the  bow  and 
stern  the  flow  is  not  two-dimensional,  as  it  is  near  amidship,   and  thus 
the  use  of  strip  theory  is  not  as  valid.     For  example,  near  the  ends  the 
water  can  vent  itself  not  only  at  the  sides,  but  also  at  the  ends,  which 
should  reduce  the  values  found  there.     These  factors  may  have  caused  the 
theoretical  values  to  be  higher  than  the  values  found  experimentally. 
As  in  heave,   it  is  believed  that  much  finer  spacing  of  sections 
should  be  used,  especially  near  the  bow  and  stern  in  order  to  obtain  bet- 
ter results  for  the  theoretical  valuesc     This  will  not  correct  the  error 
resulting  from  three  dimensional  effects,   but  would  reduce  errors  caused 
by  too  large  an  integrating  interval  in  the  vicinity  of  these  critical 
regions. 
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Chapter  VI 
CONCLUSIONS 

6.1  Evaluation  of  Step  Method 

It  is  considered  from  the  results  of  this  work  that  the  step  dis- 
turbance method  is  a  good  method  for  obtaining  added  mass  and  damping 
coefficients.  It  has  several  advantages  in  that  the  experiment  is  very 
easy  to  set  up.  By  use  of  a  device,  such  as  a  magnet,  it  requires  only 
low  potential  devices  for  power  requirements,  which  are  readily  avail- 
able. The  calculations  involved  would  be  virtually  impossible  if  done 
by  hand,  but  they  are  readily  adaptable  to  computer  operation.  One 
problem  that  may  arise  with  some  models  is  the  problem  of  wall  reflec- 
tions in  a  narrow  tank.  This  was  not  a  factor  in  the  Mariner  tests, 
but  by  orienting  the  model  to  the  centerline  of  the  tank  at  various 
angles,  this  problem  can  be  overcome. 

It  is  concluded  further  that  the  process  is  better  suited  for 
added  mass  than  for  damping.  This  is  a  factor  of  the  experiment  and 
the  analysis  program.  Damping  variations  were  much  greater  because  of 
vibrations  or  not  having  the  exact  starting  point  in  the  data*  From 
the  analysis  program,  it  was  seen  even  with  exact  generated  data,  damp- 
ing showed  variation  from  the  actual  values.  However,  if  care  is  taken 
to  get  minimum  vibrations  and  if  the  data  is  sampled  at  fine  enough  in- 
tervals, it  is  believed  that  the  process  is  valid  for  damping  as  well 
as  for  added  mass. 

A  limiting  feature  of  the  experiment  is  the  dependence  of  the  re- 
sults on  the  model's  natural  frequency.  This,  coupled  with  the 
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interference  of  the  supporting  structure's  natural  frequency,   could 
limit  the  range  of  reliable  results.     By  using  a  very  stiff  heave  rod 
and  bracing  it,   its  natural  frequency  could  be  made  high  enough  to  not 
interfere  with  the  results. 

6.2  Experimental  Results 

It  is  observed  from  the  experiments  that  both  added  mass  and 
damping  are  speed  dependent.  However,  in  the  vicinity  of  the  natural 
frequency  of  the  model,  this  dependency  is  very  slight.  The  speed  de- 
pendence was  generally  greatest  at  the  lower  frequencies. 

Virtual  mass  was  affected  by  the  model's  natural  frequency,  while 
damping  showed  little  dependence.  Thus  it  appears  possible  to  extend 
the  useful  range  of  frequencies  in  damping  by  testing  the  model  at  vari- 
ous natural  frequencies.  For  virtual  mass  calculations  the  proper 
natural  frequency  desired  must  be  used  for  the  model. 

6.3  Theoretical  Conclusions 

The  method  of  mapping  a  cylinder  into  a  ship  section  is  very  good 
if  sufficient  coefficients  are  used.  With  the  exception  of  highly  ir- 
regular sections,  almost  perfect  agreement  between  mapped  and  actual  sec- 
tions can  be  obtained.  However,  the  extent  of  the  calculations  make  it 
almost  impossible  to  perform  them  by  hand,  but  they  are  readily  adapt- 
able to  computer  operation. 

The  virtual  mass  and  damping  coefficients  derived  from  this  method 
were  believed  to  be  accurate,  as  evidenced  by  the  results  obtained  in 
heave.  However,  to  apply  strip  theory  to  the  results  in  pitch,  it  is 
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recommended  that  finer  spacing  be  used  at  the  bow  and  stern,  especially 
for  irregular  sections  such  as  are  on  the  Mariner.  By  this  method  a 
more  accurate  representation  of  the  virtual  inertia  and  damping  coef- 
ficients at  the  ends  of  the  model  could  be  obtained.  Thus,  when  multi- 
plied by  the  distance  to  the  center  of  gravity  squared,  these  sections 
would  probably  not  produce  the  effect  of  having  the  bow  and  stern  domi- 
nate in  obtaining  the  coefficients  of  virtual  inertia  and  damping  for 
the  model. 
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Chapter  VII 
RECOMMENDATIONS 

7ol  Experimental  Technique 

The  major  problem  faced  in  analyzing  the  data  was  obtaining  the 
data  in  digital  form.  The  analog  to  digital  program  used  was  very  ef- 
ficient, but  there  were  several  mistakes  made  by  the  operators  which, 
if  not  made,  would  have  greatly  aided  the  obtaining  of  the  digital 
output. 

The  first  problem  was  recording  data  at  too  slow  a  speed  and  not 
being  able  to  record  the  digitizing  signal.  The  tape  recorder  has  to 
run  at  a  speed  of  thirty  feet  per  second  to  record  a  one-kilocycle-per- 
second  square  wave  properly. 

Each  run  took  approximately  twenty- two  seconds  to  digitize.  If 
data  were  recorded  with  this  spacing  between  runs,  the  digitizing  con- 
sole should  be  set  up  to  run  continuously,  and  there  would  be  no  lost 
time  in  setting  the  start  of  each  run  when  the  computer  was  ready  for 
more  data. 

In  setting  up  the  digitizing  console,  one  has  to  set  up  an  am- 
plifying and  an  offset  section  so  as  to  keep  the  data  between  plus  and 
minus  eight  volts.  It  is  best  to  get  as  wide  a  spread  as  possible  in 
order  to  get  the  maximum  resolution  of  the  data.  In  order  to  get  this 
spread  and  not  have  any  data  be  clipped,  it  is  essential  that  all  runs 
be  recorded  to  the  same  scale.  Thus  in  setting  up  the  recording  equip- 
ment, the  zero  position  and  maximum  deflection  position  should  be  made 
the  same  on  the  recording  equipment  for  every  run0  Otherwise  reduced 
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scales  must  be  used  to  keep  all  data  within  the  allowable  range,  with 
some  loss  in  accuracy. 

The  final  point  to  be  made  about  digitizing  concerns  the  treat- 
ment of  the  magnet  parting  signal.  The  method  employed  by  the  opera- 
tors in  this  instance  produced  several  problems  that  could  have  been 
avoided.  The  first  problem  was  that  in  putting  the  one  kilocycle 
square  wave  on  the  tape,  the  instant  of  magnet  parting  was  missed  in 
some  cases  as  the  square  wave  was  put  on  the  tape  too  late.  These 
runs  had  to  be  digitized  a  second  time  to  get  useful  data.  This  was  a 
waste  of  computer  time  that  was  necessary  to  get  the  maximum  amount  of 
data  analyzed,,  Finally,  in  scanning  the  record  of  the  parting  signal, 
it  is  possible  that  the  starting  point  may  have  been  missed  by  a  couple 
thousandths  of  a  second.  This  same  accuracy  could  have  been  obtained 
by  simpler  methods.  The  best  method  appears  to  be  to  wire  the  square 
wave  generator  between  the  recorder  and  the  magnet  such  that  when  the 
magnet  parts,  the  square  wave  begins  recording.  This  process  should 
have  negligible  time  delay  „ 

The  other  big  problem  was  that  of  random  vibrations  in  the  data. 
As  stiff  a  supporting  member  as  possible  should  be  used  for  supporting 
the  ship  and  supporting  the  magnet.  This  would  not  only  reduce  the 
vibrations,  but  also  would  not  limit  results  because  of  the  equipment's 
natural  frequency.  Also,  the  model  should  be  accelerated  up  to  speed 
as  slowly  as  possible  to  minimize  the  vibrations  caused  by  this  accel- 
eration. 

To  increase  the  validity  of  the  results  found  in  heave,  a  device 
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employing  rotary  joints  and  stiff  members  may  have  the  benefit  of  less 
friction  than  the  heave  rod.  The  inertia  of  the  mechanism  may  be  dif- 
ficult to  determine,  but  it  may  be  a  better  mechanism  for  this  type  of 
experiment. 

7.2  Future  Studies 

The  method  employed  here  is  considered  a  good  one  and  could  be 
applied  to  further  studies  in  ship  motions.  It  could  possibly  be  used 
to  obtain  the  values  of  the  coupling  coefficients  between  pitch  and 
heave.  This  could  be  done  with  a  similar  arrangement  in  which  the  model 
were  allowed  both  types  of  freedom,  and  the  two  motions  recorded.  The 
method  could  also  be  extended  to  further  studies  of  damping.  This  could 
be  done  by  measuring  the  relation  between  wave  amplitude  at  the  side 
and  the  amplitude  of  the  model's  oscillation  in  heave.  The  same  linear- 
ity assumptions  could  be  employed,  and  the  results  compared  with  those 
obtained  by  forced  oscillation  tests. 
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APPENDIX  A 
DISCUSSION  OF  EQUIPMENT  FRICTION 

In  choosing  the  components  of  the  towing  apparatus  that  allowed 
pitch  and  heave  motions  of  the  ship,  and  in  choosing  the  measuring  de- 
vices, it  was  desired  to  eliminate  or  reduce  friction  to  a  minimum. 
Before  any  connection  or  measuring  plans  were  made,  this  was  the  first 
co  ns  ide  ra  tio  n. 

The  problem  in  pitch  was  not  very  great.  The  only  moving  con- 
nection, which  also  contained  the  measuring  device,  was  the  pitch 
bearing.  It  is  a  small  bearing  rotary  device  and  contributed  no  no- 
ticeable friction, 

3h  heave  the  problem  was  much  greater.  The  first  idea  was  a 
mechanism  employing  straight  bars  connected  by  small  bearing  pivots. 
While  there  would  be  several  such  pivots,  each  would  have  essentially  no 
friction,  and  the  net  effect  would  approach  frictionless  operation.  The 
biggest  problem  was  the  system  complexity  required  to  insure  heave  motion 
only.  Also,  there  was  a  problem  of  determining  the  exact  system  inertia 
in  motion. 

Because  of  the  above  problems,  it  was  decided  to  investigate  the 
use  of  a  heave  rod.  Previous  heave  rods  employed  had  a  large  amount  of 
friction.  A  new  one  was  constructed  using  a  square  rod  moving  on  a 
system  of  small  roller  bearings.  Friction  tests  were  conducted  on  it  by 
a  weight-dropping  method.  Various  weights  were  attached  to  the  rod,  and 
the  acceleration  was  measured  as  the  rod  fell.  Comparing  this  acceleration 
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to  one  gravity  unit5  one  could  determine  the  friction  in  the  mechanism,, 
Furthermore^  the  weights  were  hung  at'  fifteen  and  thirty  degrees  to  the 
vertical  to  give  the  effect  of  the  horizontal  force  induced  by  the  model 
going  through  the  water  at  various  speeds.  The  results  of  these  tests 
are  tabulated  in  Table  X, 

Table  X 
HEAVE  ROD  FRICTION 
Angle  of  Weight         Average  Friction 
0°  .385  lbs, 

15°  «506  ibs0 

30°  .822  lbs. 

To  gain  insight  as  to  the  effect  of  this  friction.,  a  test  was  con- 
ducted at  zero  speed  with  the  heave  rod  removed.  The  model  was  attached 
to  the  towing  carriage  such  that  it  was  free  to  oscilldte  on  the  linear- 
syn  only.  It  was  set  to  the  same  amplitude  as  was  used  in  the  heave 
tests0  When  the  model  was  released,,  and  the  record  of  the  heave  com- 
pared to  that  with  the  heave  rod5  there  was  no  apparent  difference  in 
the  response.  Accordingly  ^  the  effect  of  heave  rod  friction  was  neg» 
lected  in  the  final  calculations. 

As  this  test  was  conducted  at  zero  speed  and  as  the  weight  drop 
tests  showed  considerable  increase  of  friction  with  a  horizontal  com= 
ponent5  this  may  have  been  a  source  of  slight  error  at  higher  Froude 
numbers.  However^  this  effect  was  considered  secondary  and  was  neg- 
lected* 
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APPENDIX  B 
EFFECTS  OF  WALL  REFLECTIONS 

A  question  raised  early  in  the  testing  was  how  valid  the  results 
would  be  or  whether  the  effects  of  wave  reflections  in  a  narrow  tank 
would  be  such  as  to  render  the  response  meaningless.  Because  of  this 
possibility,  two  procedures  were  available  to  get  better  results  at 
zero  speed.  One  was  to  test  the  model  at  different  orientations  in  the 
tank,  and  if  this  proved  unsuccessful,  to  make  the  tests  in  a  larger 
medium,  such  as  a  swimming  pool. 

The  model  was  first  tested  parallel  to  the  tank.  The  heave  re- 
sponse obtained  is  shown  in  Figure  XVII.  The  model  was  then  tested  at 
an  orientation  perpendicular  to  the  tank,  and  the  response  is  shown  in 
Figure  XVIII,  Similar  results  were  obtained  for  pitch.  In  both  cases  it 
is  noted  that  the  model  oscillations  damped  to  zero  before  any  wall  re- 
flections affected  the  motion  of  the  model.  It  is  true  that  the  ex- 
tremely long  waves  returned  to  the  model  almost  instantly,  but  their 
effect  is  negligible,  and  it  is  the  wave  of  the  same  order  of  magnitude 
as  the  ship  length  with  which  we  are  primarily  concerned. 

It  is  interesting  to  note  that  there  is  a  considerably  longer 
zero  region  when  the  model  was  placed  across  the  tank  than  when  it  was 
oriented  parallel  to  the  tank.  This  indicates  that  wall  reflections 
must  be  considered  in  this  type  of  analysis,  but  fortunately,  for  the 
Mariner  model  the  oscillations  completely  damped  out  prior  to  wall  re- 
flections. Thus  the  results  of  the  analysis  were  essentially  the  same 
for  each  condition. 
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APPENDIX  C 
DETERMINATION  OF  COEFFICIENTS  OF  TRANSFORMATION 

The  determination  of  the  coefficients  of  transformation  of  the 
semicircle  in  one  complex  plane  to  the  ship  section  in  the  second  com- 
plex plane  proved  to  be  very  adaptable  to  computer  operation.  The  IBM 
709k   computer  at  the  Massachusetts  Institute  of  Technology  was  used 

-a 

for  this  purpose  in  this  project.  In  general,  the  mapped  cylinder  gave 
a  very  close  approximation  to  the  actual  ship  section.  The  method  used 
consisted  of  selecting  tentative  values  of  a-|_,  ao,  and  &t   by  reference 
to  Landweber  and  Macagno  (9).  The  computer  program  was  then  entered 
with  these  values  and  offsets  of  the  actual  ship  section.  An  iteration 
procedure  written  into  the  program  generated  new  values  of  the  coef- 
ficients and  offsets  which  were  compared  with  the  actual  ship  section* 
The  iterations  were  continued  until  a  reasonably  small  least  mean 
square  error  between  calculated  and  actual  offsets  was  obtained.  The 
maximum  deviation  used  in  this  project  occured  at  station  one.  Fig- 
ure XIX  in  this  appendix  compares  the  actual  station  one  with  station 
one  derived  from  the  transformation  coefficients. 
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figure  xix 

comparison  of  the  actual  mariner  station  one 
with  The  calculated  section 
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APPENDIX  D 
NEW  METHOD  FOR  DETERMINATION  OF  TRANSFORMATION  COEFFICIENTS 

After  the  main  body  of  work  on  this  project  was  completed,  a  new 
method  for  determining  the  transformation  coefficients,  a2n+l>  was  de~ 
vised.  This  method,  together  with  the  new  values  of  the  transformation 
coefficients  for  the  stations  of  the  Mariner  hullform,  are  presented  in 
this  appendix. 

Equation  (5>)  represents  the  conformal  transformation  of  the  unit 
circle  to  the  section  in  question.  In  equation  (5) 

r        ■«  -ie  (IT) 

Q  =   ipe  v   ' 


-iS 
z  =  x  +  iy  =  ire~  p  (l8) 


In  this  work  the  y  axis  points  down,  and  the  angles  9  and  3  are  measured 
from  the  vertical.  X  is  the  offset  of  the  actual  form,  and  y  is  the 
waterline.  Therefore,  when  8=3=  90°;  x  is  the  half -beam  of  the  sec- 
tion, and  y  =  0,  When  9  =  3  =  0  ,  x  =  0,  and  y  is  the  draft  of  the  sec- 
tion. 

By  substituting  equations  (18)  and  (17)  into  equation  (£)  and  sepa- 
rating the  real  and  imaginary  parts,  the  following  equations  are  obtained; 

N 

x  =  psine  +  >    (-i)n  J^r  sin(2n  +  x)e         (19) 

Z 1  x   '   2n+l 

n=0        P 
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N 
y  =  pcose  +  (-l)n+1     -|&jj     cos(2n  +  l)e  (20) 

^0  p 


It  is  then  necessary  to  solve  for  6  and  a2n+l  w^^-c^  satisfy  these  equa- 
tions. 

If  M  points  on  a  section  outline  are  available,   there  are  M  pairs 
of  equations  of  the  above  form,   or  a  total  of  2M  equations.     The  values 
of  ^n+l  ITIUS^  ^e  ^e  same    in    a^-l  °£  these  equations,  but  the  values  of  6 
will  be  unique  for  each  given  point. 

If  the  half -beam  and  draft  of  the  section  are  known,   an  additional 
equation  is  available.     As  mentioned  above,   at  9  =  90°, 

N 

i    a. 


half -beam  =  b  =  p  +      \  |££  (21 ) 

— «     P 

n=0     H 


At  e  =  o°, 


N_ 

n+1     a2n+l 


draft   =  H  =  p   +     >  (-l)irrJ-     Jg™  (22) 


n=0 


P 


2n-!-l 


These  results  lead  to  the  equation 

P(l  -  |)  +  :p-  (i  +  |)  +  ^  (1  -  |)  +  ....  -  o  (23) 

We  now  have  2M+1  equations  which  must  be  satisfied. 

In  the  new  method  of  solving  for  the  coefficients  apn+i>  ^ne  co~ 
ordinates  of  the  M  points,  as  well  as  the  half -beam  and  draft  of  the 
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section,   are  inputs  to  the  computer  program.     Three  alternatives  are  now 
available.     If  estimated  values  of  the  coefficients  a2n+]_  are  available, 
they  may  be  entered.     In  lieu  of  this,   if  estimates  of  the  angles  9  are 
available,   they  may  be  entered  as  inputs  to  the  program.     If  neither  the 
coefficients  nor  the  angles  8   are  available,  the  program  will  generate 
an  estimate  of  the  angles  9- 

Assuming  that  either  the  angles   8  are  inputs  or  that  values  of     6 
are  internally  generated,   the  program  solves  the  2M+1  equations  for  P 
and  a2n+2.  "to  give  a  least  mean  square  error  solution  between  the  gener- 
ated section  and  the  actual  ship  section.      The  values  of  the  coeffic- 
ients are  then  used  with  input  values  of  x  and  y  to  determine  a  new 
value  of  9  for  each  point.      The  program  then  re-enters  the  least  mean 
square  process  and  iterates  until  the  desired  number  of  iterations  is 
completed.     If  values  of  a^  +*    are  inputs,   the  program  computes  values 
of  ©  before  entering  the  mean  square  process. 

It  should  be  noted  that  the  program  requires  the  number  of  co- 
efficients desired  and  the  number  of  iterations  desired  as  inputs. 

As  stated  above,  the  program  finds  values  of  p.  However,  since 
we  are  interested  in  the  case  of  the  unit  circle,  the  program  is  nor- 
malized with  p  equal  to  one. 

Using  this  new  method,   new  values  of  ^n*-!  were  determined  for  the 
Mariner  hullform.     These  coefficients  are  presented  in  Tables  XI  and  XII. 

It  is  believed  that  this  improved  approximation  to  the  hullform 
will  not  appreciably  affect  the  values  of  damping  and  virtual  mass  for 
the  Mariner.     However,   in  future  work  it  is  recommended  that  the  new 
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method  presented  above  be  used. 

A  listing  of  the  analysis  program  in  MAD  computer  language  is  at- 
tached to  this  appendix. 


Table 

XI 

TRANSFORMATION 

COEFFICIENTS 

Station 

!i 

!2 

!£ 

!l 

ao_ 

0 

-.858571 

-.118183 

-.018627 

+.006032 

-.000315 

1 

-.675965 

-.085363 

-.513531* 

+.006935 

+.005787 

2 

-.h79528 

-.0U3086 

+  .0101*33 

+.003225 

+.013838 

3 

-.286588 

-.022750 

+  .0191*79 

+.003222 

+.016652 

a 

-.117882 

-.020811 

+.019821 

+.005278 

+.010582 

S 

-.000^35 

-.035251 

♦.017899 

+.002628 

+  .00781*1* 

6 

+.076025 

-.053853 

+.016527 

+.001157 

+.005735 

7 

+.121668 

-.073578 

+  .01071*7 

-.000671* 

+.002915 

8 

+  .ll*399l* 

-.102216 

♦  .  001*109 

-.001911 

9 

+.150007 

-.1190U0 

-.001388 

-.002853 

10 

+.150031 

-,12li7U6 

-.OOU275 

-.001*280 

+.001686 

11 

+  ,  151110* 

-.122959 

-.0031*35 

-.ooiaoii 

12 

+.150825 

-.110575 

-.001117 

-.001822 

13 

+.1U8566 

-.085017 

+.005537 

-.001*659 

Hi 

+.139671 

-.OU8057 

+.016075 

-.001*566 

+.001582 

15 

+.1171*10 

-.0021*30 

+  .02)431*5 

-.001183 

+.006575 

16 

+.O69663 

+,053825 

+.032166 

+  .00001*2 

+o00933U 

17 

-.02691*9 

+.106876 

+.0h35l2 

-.001900 

+.006618 

18 

-.216518 

+  .ll;6021 

+  .071*115 

+.012576 

-.003023 

19 

-.562587 

+.131206 

+  .09821*5 

+cOU7390 

+.016198 

20 

+.227269 

+.251236 

-.121012 

+.OU3271 

-.109952 

A-ll 

i 

rable  XII 

COMPARISON  OF  ERRORS 

WITH  THOSE 

FOUND  BY  PROCESS 

USED  IN  THE  MAIN  BODY  OF 

THE  PROJECT 

New  Method 

Main 

Body- 

Station 

Average 
Excursion 

Maximum 
Excursion 

Average 
Excursion 

Maximum 
Excursion 

0 

0.07 

0.11 

1 

0.09 

0.17 

0.32 

0.93 

2 

0.16 

0.33 

0.17 

0.32 

3 

0.13 

0.25 

0.20 

o.$k 

h 

0.12 

0.23 

0.11 

o.ho 

5 

0.10 

0.22 

0.13 

0.32 

6 

0.08 

0.17 

0.08 

0.17 

7 

0.06 

0.13 

0.06 

0.18 

8 

0.03 

0.07 

0.07 

0.13 

9 

0,02 

0.07 

0.02 

0.08 

10 

0.05 

0.12 

0.08 

0.23 

11 

0.06 

0.11 

0.06 

0.11 

12 

0.03 

0.09 

0.03 

0.09 

13 

0.03 

0.05 

o.oU 

0.08 

1U 

0.03 

0.09 

o.oU 

0.10 

15 

0.09 

0.25 

0.19 

o.Uo 

16 

0.08 

0.19 

0.2li 

o.55 

17 

0.12 

0.27 

0.19 

0.52 

18 

0.06 

0.11 

0.22 

0.U6 

19 

0.15 

0.30 

0.29 

0,69 

20 

0.07 

0,13 

0.07 

0.20 
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RCOEFFICIENT  DETERMINATION  FOR  HULL  FORMS 
DIMENSION  B(625.BDIM)  ,THETA(625,TDlM)»X(30)»Y(30)iA6AR(2t>)» 

1  APRI(15)»BETA(30)tCRlT(20)»A(15)»WYE(30).tX(30).LRRX(30)t 

2  ERRY(30  )  ,NAME(  12)  » DEGREE  (  62  i> .  00  I  M  )  »ALPHA(  30)  »DUM<  ^2^»DIMD) 
INTEGER  I  ,  J »K,M,N»Q»COUNT,PWR, SIGNAL 

VECTOR  VALUES  BDIM  =  2,1,25 
VECTOR  VALUES  TDIM  =  2,1,25 
VECTOR  VALUES  DDIM  =  2,1,25 
VECTOR  VALUES  DIMD  =  2,1,15 
START      READ  FORMAT  T I TL1  ,NAME ( 1 ) • . . NAME  I  1 2 ) 
VECTOR  VALUES  TITL1  =  $(12C6)»$ 
PRINT  FORMAT  T I TL2 , NAME ( 1 ) . . .NAME ( 1 2 ) 
VECTOR  VALUES  TITL2  =  $ t 1H1 » 12C6 ) *$ 
READ  FORMAT  INPUT,  H ,m ,n ,0 , BEAM , S I GN AL 
VECTOR  VALUES  INPUT  =  i ( F10. 6 , 3  I  2 , F 10. 6 , I  2 ) *$ 
BDIM(2)  =  N+2 
DIMD(2)  =  N+l 
TDIM(2)  =  0+1 
DDIM(2)  =  0+1 

READ  FORMAT  OFFSET,  X(l).#.XtM) 
VECTOR  VALUES  OFFSET  =  $(5F10.4)*$ 
READ  FORMAT  WL ,  Y(1)...Y(M) 
VECTOR  VALUES  WL  =  $(5F10.4)»$ 
KAY  =  BEAM/H 
El  =  1.0E-06 
WHENEVER  SIGNAL  .E.  0 
READ  FORMAT  ANGLE,  EX ( 1 ) . . • EX ( M) 
THROUGH  BEE1,  FORM  =  1,1,  I.G.M 
BEE1       DEGREE( 1,1 )  =  EX( I ) 

VECTOR  VALUES  ANGLE  =  $(5F10.5)*$ 

THROUGH  TWENTY,  FOR  I  =  ltltI«G«M 
TWENTY      THETA(Itl)  =  DEGREE ( I , 1 ) /57 .295780 

TRANSFER  TO  TWELVE 
OR  WHENEVER  SIGNAL  .E.  1 
READ  FORMAT  AV/L,  AH)...AIN) 
VECTOR  VALUES  AVAL  =  $ ( S4,4E14.6 ) *$ 

THROUGH  THIRTY,  FOR  I  =  1 , 2  ,  I .G. ( 2*N-1 ) 
THIRTY     A(2»I-1 )  =  A( I ) 
DEN  =  1.  +  A(l ) 

THROUGH  FORTYt  FOR  I  =  3  ,2 , I .G. ( 2*N-1 ) 
FORTY       DEN  =  DEN  +  At  I ) 

RO  =  BEAM/DEN 

THROUGH  FIFTY,  FOR  I  =  1  ,2 ♦ I .G. ( 2*N- 1 ) 
FIFTY      APRI ( I )  =  At  I )*R0 

TRANSFER  TO  TWFLVF 
OTHERWISE 

CONTINUE 
END  OF  CONDITIONAL 
R 
R 
THROUGH  INITL,  FOR  I  =  1,1, I.G.M 
THETA(Itl)  =  ATAN. (X( I )/Y( I ) ) 
INITL      WHENEVER  Y(I)  .LE.  El,  THETA(Itl)  =  1.57080 
TWELVE     THROUGH  ONE,  FOR  J  =  1,1,J.G.Q 

WHENEVER  SIGNAL  .E.  1  .AND.  J  .E.  1,  TRANfFER  TO   Bl 


\ 


TWO 
THREE 


THROUGH  TWO.  FOR  I  =  1.1.  I.G.M 

R  (  I  .1  )  =  SIN.  ( THETA(  I  ,j)  ) 

B(  I  .2)  =  B(  I  .1  ) 

B( I  » (N  +  2  )  )  =  X(  I  ) 
B((I+M).l)  =  COS.  (THETA(  I  tJ)  ) 
B(  (  I+M)  .2)  =  -F (  (  I+M)  ,1 ) 
B(  (  I+M)  ♦ (N+2)  )  =  Y(  I  ) 

THROUGH  TWO.  FOR  K  =  3.1.  K.G.(N+1) 
B(ItK)  =  (  (-1 . ) .P.K  )*SIN. (  (  ( 2*K)-3 )*THETAI  I »J )  ) 
B((I+M).K)  =  -(  (-1. ).P.K)*COS.(  (  (2*K)-3)*THETA(  I »J)  ) 
THROUGH  THREE*  FOR  K  =  1 . 1 . K .G. ( N+l ) 
B((2»M+1).K)  =  1.  +  ( (-1, ).P.K)*KAY 
B( (2*M+1) »(N+2) )  =  0, 


EXECUTE  GLSOH. ( B.ABAR,DUM. (N  +  l  ) ♦ (2*M+1  )  ) 
PRINT  RESULTS  ABAR ( 1 ) . . . ABAR ( N+l ) 


FOUR 


FIVE 


SIXPRI 
SIX 


EIGHT 


RO  =  ABAR(1 ) 

THROUGH  FOUR,  FOR  I  =  2.1. 
APRI  (  (2*1  )-3)  =  ABAR( I ) 


I.G.  (N+l) 


THROUGH  FIVE.  FOR 
A( I )  =  APRI ( I ) /R6 


I  =  1.2.  I .G.(  (2»N)-1  ) 


FOR  I  =  1.1.  I.G.M 

+  APRI  ( 1 )  )»SIN. ( THETA(  I  .J)  ) 


THROUGH  SIX. 

EX(I)  =  (RO 

PWR  =  0 
DEGREE(I.J)  =  57.2958*THETA( I .J) 

WYE(I)  =  (RO-APRI ( 1 ) )*COS. ( THETA( I .J) ) 

THROUGH  SIXPRI.  FOR  K  =  3.2.  K.G. < ( 2*N ) -1 ) 

PWR  =  PWR  +  1 
WYE(I)  =  WYE(I)  -  ( (-l.).P.PWR)*APRI (K)*COS. (K*THETA( I.J) ) 
EX(I)  =  EX(I)  +  ( (-1. ) .P. PWR )*APRI (K)*SIN. (K*THETA ( I . J) ) 
FRRY(I)  =  .APS.(WYE(I)  -  Y(I)) 

ERRX(  I  )  =  .ABS. ( EX ( I  )-X (  I )  ) 


ALPHA(O)  =  0, 

BFTA(O)  =  0. 

THROUGH  EIGHT.  FOR  I  =  1.1.  I.G.M 

ALPHA(I)  =  ALPh'A(I-l)  +  ERRX  (  I  )  *ERRX  (  I  ) 

BETA(I)  =  BETA(I-l)  +  ERRY ( I  )  *ERRY (  I  ) 


CRIT(J)  =  SORT. ( ALPHA(M)  +  BETA(M)) 

EXECUTE  OUTPT, (J.RO.A.CRlT.X.EX.ERRX.Y.WYE.ERRY.DEGREE.M.N) 
WHENEVER  J.E.Q.  TRANSFER  TO  FINIS 


THROUGH  CEE1,  FOR  I  =  lt2t  I.G.(2*N-1) 
CFE1         PRINT  FORMAT  ARC »  A{T) 

VECTOR  VALUES  ABC  =  $(1H  »E14.6)*$ 
PRINT  FORMAT  LMS»  CRIT(J) 

VECTOR  VALUES  LMS  =  $(21H  LEAST  MEAN  SQUARE  =  .F14.6)*$ 
PRINT  COMMENT  SO  MEASURED  Xd)        CALCULATED  X(I)        X(I) 
1  ERROR         MEASURED  Y(I)        CALCULATED  Yd)        Yd)  ERR 
20R       THETA(I)$ 
THROUGH  TEN.  FOR  I  =  Ltlt  I.G.M 
TEN  PRINT  FORMAT  OFVAL*  X ( I ) , EX ( I ) , ERRX d ) » Y ( I ) ♦ WYE ( I ) . ERRY ( I ) , 

1   DEGREE(I.J) 

VECTOR  VALUES  OFVAL  =  $(1H  t F 14. 6 »S5 » F 1 4 . 6 . S5 » F14 . 6 » S5 t 
1   F14.6tS5tF14.6tS5tF14.6tS5tFl4.6)*$ 

FUNCTION  RETURN 

END  OF  FUNCTION 

MAD 

EXTERNAL  FUNCTION  (A.XtBtNtM) 

ENTRY  TO  GLSOH. 

INTEGER  ItJtKtNtM 

THROUGH  QQlt  FOR  I  =  ltl»  I.G.N 

X( I )  =  0. 

THROUGH  QQlt  FOR  J  =  ltlt  J.G.M 
001        X(I)  =  X(I)  +  A(  Jtl  )*A(  Jt(N+l)  ) 

THROUGH  QQ2t  FOR  I  =  1»1»  I.G.N 

THROUGH  0Q2t  FOR  J  =' 1.1»  J.G.N 

B(  I  »J)  =  0. 

THROUGH  Q02»  FOR  K  =  ltlt  K.G.M 
QQ2         B(ItJ)  =  BdtJ)  -+  A(K.I }*A(K»J) 

THROUGH  QQ3.  FOR  I  =  1»1»  I.G.N 

WHENEVER  B( I . I )  .E.  0. 

THROUGH  QQ4»  FOR  K  =  (I+l)tl«  K.G.N 

THROUGH  QQ5»  FOR  J  =  ltlt  J.G.N 

TEST  =  B( I tJ) 

B(  I  tJ)  =  B(KtJ) 
Q05        B(KtJ)  =  TFST 

TEST  =  X( I ) 

X(I  )  =  X(K) 

X(K)  =  TEST 
QQ4        WHENEVER  B(Itl)  .NE.  O.t  TRANSFER  TO  QQ6 

END  OF  CONDITIONAL 
0Q6        STORE  =  B( I t I ) 

THROUGH  QQ7t  FOR  J  =  ltlt  J.G.N 
QQ7        B( I tJ)  =  B( I tJ) /STORE 

X(  I  )  =  X(  I  )/STORE 

THROUGH  QQ3t  FOR  K    =  ltlt  K.G.N 

WHENEVER  K.E.I.OR.BKtl  )  .E.O.t  TRANSFER  TO  QQ3 

STORE  =  B(KtI ) 

THROUGH    QQ8t     FOR    J    =     Itl.t     J.G.N 

B(KtJ)     =    BdtJ)*STORE    -    B(KtJ) 
QQ8  WHENEVER    K.L.It     B(KtJ)     =    -B(KtJ) 

X(K)     =    X(  I  )*STORE    -    X(K) 

WHENEVER  K.L.It  X(K)  =  -X(K) 
QQ3        CONTINUE 

FUNCTION  RETURN 

END  OF  FUNCTION 
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